OU-HET 555 
WU-AP/240/06 
hep-th/0602242 



Inflation from Superstring/M-Theory Compact iflcat ion with 

Higher Order Corrections II 

— Case of Quartic Weyl Terms — 

Kenta Akune,*^^ Kei-ichi Maeda"'*'"^'^ and Nobuyoshi Ohta"^'^ 



Department of Physics, Waseda University, Shinjuku, Tokyo 169-8555, Japan 
^Advanced Research Institute for Science and Engineering, Waseda University, Shinjuku, Tokyo 169-8555, 

Japan 

Waseda Institute for Astrophysics, Waseda University, Shinjuku, Tokyo 169-8555, Japan 
'^Department of Physics, Osaka University, Toyonaka, Osaka 560-0043, Japan 



Abstract 

We present a detailed study of inflationary solutions in M-theory with higher order quantum corrections. 
We first exhaust all exact and asymptotic solutions of exponential and power-law expansions in this theory 
with quartic curvature corrections, and then perform a linear perturbation analysis around fixed points for 
the exact solutions in order to see which solutions are more generic and give interesting cosmological models. 
We find an interesting solution in which the external space expands exponentially and the internal space is 
static both in the original and Einstein frames. Furthermore, we perform a numerical calculation around this 
solution and find numerical solutions which give enough e- foldings. We also briefly summarize similar solutions 
in type II superstrings. 



^e-mail address: akune@gravity.phys.waseda.ac.jp 
^e-mail address: maeda@waseda.jp 

^e-mail address: ohta@phys.sci.osaka-u.ac.jp. Address after 31 March 2006: Department of Physics, Kinki University, Higashi- 
Osaka, Osaka 577-8502, Japan 



1 Introduction 



The recent cosmological observations have confirmed the existence of the early inflationary epoch and the 
accelerated expansion of the present universe [1]. Though it is not difficult to construct cosmological models 
with these features if one introduces scalar fields with suitable potentials, it is desirable to derive such a 
model from fundamental theories of particle physics without making special assumptions. The most promising 
candidates for such theories are the ten-dimensional superstrings or eleven-dimensional M-theory, which are 
hoped to give models of accelerated expansion of the universe; upon compactification to four dimensions. 

There is a no-go theorem which forbids such solutions if six- or seven-dimensional internal space is a time- 
independent nonsingular compact manifold without boundary [2] . This theorem also assumes that the gravity 
action does not contain higher curvature corrections. So we can evade this theorem by violating some of the 
assumptions. In fact, it has been shown that a model with certain period of accelerated expansion can be 
obtained from the higher-dimensional vacuum Einstein equation if one assume a time-dependent hyperbolic 
internal space [3] . It has been shown [4] that this class of models is obtained from what are known as S-branes 
[5, 6, 7] in the limit of vanishing flux of three- form fields (see also [8]). For other attempts for inflation in the 
context of string theories, see, for instance, Refs. [9, 10, 11]. 

The scale when the acceleration occurs in this type of models is basically governed by the Planck scale 
in the higher ten or eleven dimensions. With phenomena at such high energy, it is expected that we cannot 
ignore quantum corrections such as higher derivative terms in the theories at least in the early universe. It 
is known that there are terms of higher orders in the curvature to the lowest effective supcrgravity action 
coming from superstrings or M-theory [12, 13, 14, 15, 16]. The no-go theorem does not apply to theories with 
higher derivatives because they violate the assumptions. With such corrections, they will significantly affect 
the inflation at the early stage of the evolution of our universe. 

The cosmological models in higher dimensions were studied intensively in the 80's by many authors [17, 
18, 19, 20, 21]. It was shown that inflation is indeed possible with higher-order curvature corrections [19, 20]. 
However, most of the work considered theories with higher orders of scalar curvature, which are not typical 
correction terms known to arise in the superstring theories or M-theory. It is thus important to examine if 
the above result of small e-folding is modified with higher-order corrections expected in these fundamental 
theories. 

In particular, the leading quadratic correction for heterotic string theories is proportional to the Gauss- 
Bonnet (GB) combination [12, 13, 14]. This model has been studied in some detail in Ref. [19] and it was shown 
that there are two exponentially expanding solutions, which may be called generalized de Sitter solutions since 
the size of the internal space also depends on time. Despite the exponential behavior in the original frame, 
those solutions give non-inflationary power-law expansion in the Einstein frame as indicated in our previous 
papers [22, 23]. On the other hand, there is no quadratic and cubic curvature corrections in type II superstring 
theory or M-theory, and thus the first higher order corrections start with quartic curvature terms [15, 16]. 

In our previous papers [22, 23], we have reported our results on this problem for M-theory with forth- 
order terms in terms of the Riemann curvature tensors. However, it turns out that the coefficients we took 
had opposite signs to the M-theory case. Since there is ambiguity in the coefficients coming from the field 
redefinition, this does not immediately invalidate the results, but it is necessary to study how the results 
change. Here we give the details of our results in M-theory as well as type II superstrings with the correct 
coefficients. We consider these corrections given in terms of the Weyl tensors which are favorable because 
only corrections in this scheme do not affect the highly symmetric tree-level solutions such as AdSy x 5^ [16]. 
We also inchide an additional quartic term in scalar curvature with an arbitrary coefficient S in order to take 
into account the ambiguity mentioned above. With this action, we exhaust exact solutions as well as past 
and future asymptotic solutions and then discuss inflationary solutions among them. The past and future 
asymptotic solutions are useful in describing the inflation at the early universe and the present accelerating 
cosmology, respectively. Furthermore, we perform a linear perturbation analysis in order to see which solutions 
are more generic and to make interesting cosmological models. We find an interesting solution in which the 
external space expands exponentially and the internal space is static both in the original and Einstein frames. 
This may be regarded as "moduli stabilization" by higher order corrections, but this is not the usual moduli 
stabilization in the sense that the moduli are fixed and stable. Rather we are interested in such solutions in 
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which the sizes of the internal spaces do not grow too much while they exhibit inflation. What we find is that 
this is possible with higher order corrections. We also perform a numerical calculation around this solution 
and find that some spacetimes give enough e- foldings. Finally we briefly summarize similar solutions in type 
II superstrings. Necessity of the higher order corrections for inflation in M-theory is also discussed in Ref. [24]. 

In the next section, we present our actions and field equations to be solved. We write down these for 
D — {1 + P + q) dimensions with p external and q internal space dimensions. Though we are mainly interested 
in p = 3 in this paper, there may be interesting applications if we keep the dimension p arbitrary. We give the 
equations for maximally symmetric spaces with non-vanishing curvatures. The explicit forms of the actions 
are given in Appendix A, and the field equations in Appendices B and C. Although similar equations are given 
in our previous papers [22, 23], the present theory is different in that the quartic terms are written in the Weyl 
tensor instead of Riemann, and also that we have an additional term R^. So we present these equations for 
our new system. 

We give exact solutions as well as past and future asymptotic solutions with exponential and power-law 
expansions in § 3 for 5 = and in § 4 for ^ ^ 0. Also we present the solutions for maximally symmetric spaces 
with non- vanishing curvatures. 

In § 5, we perform a linear perturbation analysis around fixed points for the exact generalized de Sitter 
solutions that are given in § 3 and § 4 in order to see which solutions are more generic and to make interesting 
cosmological scenario. 

In § 6, we also summarize exact solutions for type IIB superstrings for constant dilaton. 

Using the obtained solutions, we discuss an inflationary scenario in §7. Many of our exact solutions do not 
seem to give successful inflation in the sense that they do not give big enough e- foldings. However the simple 
analysis of exact solutions does not tell us what happens after the inflationary solutions decay. Actually it 
turns out by numerical analysis that there are interesting solutions, which first approach to the exact solution 
for which the external space expands exponentially and the internal space is static both in the original and 
Einstein frames, and then eventually go to a stable solution. For such solutions, we show that it is possible to 
obtain enough e-foldings for successful inflation. This is a very interesting possibility of achieving inflationary 
solutions. 

The contents of §§ 3 - 6 summarize our analysis of solutions and are rather technical. The reader may 
find it useful to skip this part for the first reading, and then come back to check when they study physically 
interesting solutions described in § 7. 

2 Field equations 

We consider the low-energy effective action for M-theory {D = 11) with higher order corrections keeping 
dimension D arbitrary: 
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(2.1) 



with 




D 



(2.2) 




(2.3) 
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(2.5) 
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Sr^ = ^Jd''x^R\ (2.7) 

where 

^2n — 2"(£' — 2n)! tQ,j...Q.Q_2„pi<Ti...p„cr„-'i- /iii^i n„i'n t V^-°/' 

and i?* is a quartic term of scalar curvature. Here we have dropped contributions from form fields, is a 
_D-dimcnsional gravitational constant, and ai, . . . , 0:4, 7 and 5 axe numerical coefficients. The coefiicient ai of 
the Einstein-Hilbert (EH) term is 1 by definition, and though is zero for all superstrings and M-theory, we 
have included it since it will be useful for examining other cases. The Weyl tensors in L^/ are defined by 

R 

{gxvQuK - 9Xk9hv)- (2.10) 



{D-l){D-2) 



In our previous papers [22, 23], we gave field equations for this system with arbitrary couplings, but we 
focused on the case of nonzero and 7. We discussed the flat external and internal spaces in Ref. [22] with 
various coefficients of the higher order terms, whereas all combinations of curved spaces are examined in [23]. 
We considered the case in which the quartic term Lw are written in terms of Riemann tensors. However, it 
turned out that the coefiicients we took in [23] were opposite in sign to the M-theory case,^ and we should set 

a, = as = 0, "4=3,/29x(2^)4' ^ = 3 x (L)^ ' ^^'^^^ 

where T2 = {2Tr^ / K^y^^ is the membrane tension. Though the field equations remain valid, the numerical 
results on the generalized de Sitter solutions are significantly affected by this sign change. We find that many of 
the solutions found in our previous paper go away if we simply reverse the signs of these coefficients. However, 
we should also note that contributions of the Ricci tensor R^^ and scalar curvature R are not included in the 
fourth-order corrections (2.9) because these terms are not uniquely fixed. This means that there is significant 
ambiguity in the additional terms involving these tensors, and in particular this allows us to put the forth- 
order terms in terms of the Weyl tensors as given above. This form appears particularly favorable because only 
corrections in this scheme do not affect the highly symmetric tree-level solutions such as AdSr x S'^ (AdSs x 
for type IIB superstring theory) [16]. In view of this situation, it appears more appropriate to consider the 
quartic correction terms given by Weyl tensors. 

Furthermore it is interesting to examine how the ambiguity may affect the results. For this purpose, we 
also include additional quartic Ricci scalar term (2.7). Here we discuss what value of 5 is plausible. Writing 
down Eq. (2.9) in terms of Riemann curvature tensor, we have the following equation: 

Lw[R^lupa) - Lw[<^iJ.upa) + (£)_ 1)2(£)_3)3-" ^"Z^"" + (D - 1)2 (£) - 3)3 

-|- (terms containing both C^^^o- and Ra0,R) , (2-12) 
where we omit terms which contain both C^^pa and Rap, R, and LwiR^upa) is defined by 

T ( T) \ tjXlivk d d puot D0 I ^ dXkuu d d pea p/3 to i q\ 

^Wy-tlfivpa) — n, ^anvl3^X ptrK + 2 ^afinv^X paK ■ {^■^'^) 

We thus find that the difference between Lw in terms of the Riemann and Weyl tensors depends on the Ricci 
and scalar curvatures with very suppressed coeflacient due to large D = 11. It is then natural to consider the 

^ See the recently replaced version of Ref. [16] . 
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i?^ correction terms with the same order of magnitude as above. Otherwise the scalar curvature terms will be 
dominant and the stringy effects may not be seen. Thus, the appropriate value of 6 appears to have the order 
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(£)- l)2(£)-3)3 



7 ~ 10 ^7 . 



(2.14) 



However, since we do not know the exact contribution from Ricci and scalar curvatures, we leave 6 to be 
free. 

This is the system that we are going to examine. 



2.1 Basic equations for cosmology 

The metric of our Z)-dimensional space is 

dsjj = -N'^{t)dt'^ + a^{t)dsl + 6^(i)rfs^ , 

with 

N{t) = e"''^*) , a{t) = e"i(*) , b{t) = e"^(') 



(2.15) 



(2.16) 



where D = 1 + p + q. The external p- and internal g-dimensional spaces {ds^ and ds^) are chosen to be 
maximally symmetric. The curvature constants of those spaces are defined by ap and Oq. The sign of (jp 
{(Tq) determines the type of maximally symmetric spaces, i.e. Up (or Oq) = —1, and 1 denote a hyperbolic 
space, a flat Euclidean space, and a sphere, respectively. The hyperbolic and flat spaces are supposed to be 
compactified by identifying boundaries of those finite part. 

From the variation of the total action (2.1), whose explicit forms with the metric (2.15) are given in 
Appendix A, with respect to uq,ui and W2, we find three basic field equations: 



F(p) 



E /i^^ + ^ E sir' +Yt: /.(f) + F^^ + F^J = , 

n=l n=l n=l 

E fi'' + ^ E 5^"' + ^ E h^r?' + + = ' 



where X = ili — uoU\ + u\, Y = U2 — uoU2 + u^, and 



Fn 


= Fn{uo,Ul,U2,Ap,Aq) , 








Fw 


= Fw{uo,ui, 


U2,Uo,Ul,U2, 


ill, 


il2, 'ill, 


U2,X,Y,X,Y) , 


Fr^ 


= Ffii{uo,ui, 


U2,Uo,Ui,U2, 


ill, 


U2, Ul 


U2,X,Y,X,Y) , 


f(p) 

J n 




U2,Ap,Aq) , 




= 

Hn 


g^P\uo,Ui,U2,Ap,Aq) , 


"ri 




U2,Ap,Aq) , 








pip) 


= fI^\uo,ui 


,U2,Uo,Ui,U2 


ill 


ll2, Ul 


,U2,X,Y,X,Y,X,Y) , 


pip) 


= fI^^{uo,ui 


,U2,Uo,Ui,U2 


ill 


,il2, Ul 


,U2,X,Y,X,Y,X,Y) , 


fil) 
J n 




U2,Ap,Aq) , 




9^' = 


9nH'^0,Ul,U2,Ap,Aq) , 


"-n 




U2,Ap,Aq) , 








pil) 


= fI^\uo,ui 


U2,Uo,Ul,U2 


ill 


il2, Ul 


,U2,X,Y,X,Y,X,Y) , 


pig) 


= fI^}{uo,Ui 


,U2,Uo,Ui,U2 


ill 


1*2, Ul 


,U2,X,Y,X,Y,X,Y) , 



(2.17) 
(2.18) 
(2.19) 



(2.20) 
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are explicitly given in Appendix B. Here Ap and Aq are defined by 

Ap = ul + cTp exp[2(Mo - Ml)], 

Ag=u2 + (Jqexp[2{uo-U2)]. (2.21) 

Since uo is a gauge freedom of time coordinate, we have three equations for two variables ui and U2. It 
looks like an over-determinant system. However, these three equations are not independent. In fact, we can 
derive the following equation after bothersome calculation: 

F + {pill + qu2 - uo)F = p-iiiF^P^ + qu^F'^'i^ . (2.22) 

The constraint equation = is satisfied if other dynamical equations are solved and it is initially satisfied. 
As argued in Ref. [23], it is in general enough to solve the two equations F = Q and F'^p^ = (or F^^^ = 0) 
instead of trying to solve all three equations. 

2.2 Ansatz for solutions 

We now analyze our basic Eqs. (2.17) - (2.19) and look for inflationary solutions. Since we are interested in 
analytic solutions, we study the following two cases: 

(1) Generalized de Sitter solutions: 

Using a cosmic time, i.e. uq = 0, an exponential expansion of each scale factor is given by Ui = /Uf + lnao, 
and U2 = vt + ln6o, where /x, v, ao and bo are constants. 

(2) Power-law solutions: 

To find a power-law solution, although we can discuss it with the above cosmic time, we use a different 
time gauge, which is defined by uq = t. Using this time coordinate, a power-law solution is given by 
ui = fit + Inao, and U2 = i^t + ln6o, where /U and u are constants. 

If we choose the following time coordinate uq, 

Uo = et, ui= iJ,t + Inao, and U2 = i^t + ln6o, (2.23) 

we can discuss both solutions in the same set up, that is, e = for case (1), while e = 1 for case (2). In the 
latter case, in terms of a cosmic time r = e*, we see that the solution gives the power-law behavior: 

a = e"i=aor'', and b = e""^ = bor" . (2.24) 

Note that when the curvature constant ap (or u,) vanishes, ao and bo are arbitrary but we can set ao = 1 (or 

6o = 1) because such a numerical constant can be absorbed by rcscaling of the spatial coordinates. 
Before giving the solution, we note on the unit used in our solutions. Rcscaling a^, 7, /x and v as 

a4 = a4/|7|, 7 = 7/l7l(=l), fl = &nd u = (2.25) 

we can always set 7 to +1. We also have to rescale time coordinate ast = \^\~^^^t. The typical dynamical time 

scale is then given by l^yl^^^ ~ 0(to^^), where mo = k~ij^^^^ is the fimdamental Planck scale. In particular, 
for M-theory, we find = 6-i/6(47r)-^/9m]"i^ ~ OASlSUGrriii from Eq. (2.11). After this scaling for the 

M-theory, we have 

"*=3^' ^ = ^^-^^^ 

and a free parameter S. 

We use the above unit throughout this paper and omit tilde for simplicity. We now present solutions for 
(5 = and J 7^ in order. 
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3 Solutions in M theory for ^ = 



Prom Eq. (C.2), we expect there may exist no exact solutions expect for the case (Jp = CTq = 0. However, 
even for the case of CTp 7^ or ^ 0, we may have some asymptotic solutions either in the future direction 
[t 00) or in the past direction {t —00), which describe the universes in these stages. We classify solutions 
for Eqs. (2.17) - (2.19) by the signatures of and cr, in the following subsections. 

3.1 Up — <Jq — Q 

In this case, Ap = ii^,Aq = are constants. We will discuss the cases of e = and e = 1 in order. 

3.1.1 Generalized de Sitter solutions (e = 0) 

In this case, we can take the independent equations as F = and F^^^ + F^'^^ = for /x, 0. From the 
explicit forms of field equations given in Appendix C, we have two algebraic equations for p = 3: 

o.\ [6/L(^ + 6(7/ii^ + q\v^\ + a^v^ [336(74/v,'' + 168g5/i^;/ + 24(7g/iz/^ + q'jv\ 

+ (/+2)4(g + 3)3 [~^(^ ~ ^^^1^*' ^^-"^ + ~ + ~ ^^^^ + + ~ ~ 1512)/x3z. 

- 2(g - l)(14g^ + 13g^ - 242g^ + 1079g^ + 4014^^ - 2970g - 12m)p?v^ 
+ 4(2g^ - 2,q^ - llgS + 413g* + 337g^ - 1494^^ + 900g + 2\<a)iiv^ + (g - 7)7Vi(3, q)v^\ 

+ ^(-12/^2 + %qiiv + q{q - 7^) [12^^ ^ q^^^ ^ ^ i)^^2j 3 ^ ^ ^ 
(/X - i^){ai + 4^4 [30{q - l)iti''u^ + 12(g - 1)5^2^' + {q- 1)6Z^^] 

- 3(g - 1)(2/ + 17g^ - llq^ - 142g^ + 1109g^ + 1995g^ - 3978g - 2376)/i^i/ 
+ {q- l)(7g^ + 40q^ - lOAq^ - lUq^ + 4887g^ + 6372^^ - 14580g - 7776)/xV2 

- 2g(/ + Aq^ - llq^ + 135/ + 1104g^ + 9g^ - 3366g + 3024)/iz/^ + 67Vi(3, q)u^] 

+ 46[l2n^ + 6qi^u+{q + l)oiy^f}=0, (3.2) 

where Ni {q, p) is an integer constant defined by 

Ni{q,p) = p\ip\q - 2) + p{q^ + lq- 14) - 2q{2q^ _ 5g + 1) - 7) 

+ {q- lf{p\q - 3) + ipq{q + 1) - Sg^) . (3.3) 

Using the values for 0:4 and 7 in Eq. (2.26) and setting q = 7, we can solve these equations numerically and 
found the following solution for ^ = 0: 

(^,1/) = ME1±(±0.10465, tO.93666) . (3.4) 

Here ME1± (the first exact solutions in M theory) are the names given to the solutions. We will use similar 
names for solutions in what follows. 

3.1.2 Power-law solutions (e = 1) 

Setting e = 1 in Eqs. (C.4) - (C.29), there is no exact solutions. We find that the EH action is dominant 
as t ^ 00, while the actions S4, Sw and 3^4. become dominant as t ^ —00. Here we present asymptotic 
power-law solutions for each case. 
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As i ^ oo, EH term dominates and our basie equations reduce to 

Pif/ + qii^"^ + 2pqfj.v = 0, (3.5) 
q:^{:,~^i-l)-{p-l)tx = 0, (3.6) 
p^ip - 1^ - 1) - {q - 1)1^ = 0. (3.7) 

We can easily show that these three equations are equivalent to the following two equations, if it is not 
Minkowski space {fi = u = 0): 

pfi^ + qv^ = 1, p/j, + qiJ = 1, (3.8) 
which is a special case of Kasner solutions. We have a solution 

P± Vm(p + 9- 1) 



p{p + q) 

qT \/pq{p + q- 1) 



(i{p + q) 

For p = 3,q = 7, we get two future asymptotic solutions: 



(3.9) 



As f ^ — oo, the fourth-order terms dominate. So let tis present asymptotic power-law solutions only with 
quartic terms. Assuming the metric (C.l) with e — I, our basic Eqs. (2.17) and (2.19) give three algebraic 
equations. For p = 3, q = 7,we have 



■ [5913(Ai - 7)(/i - l)-"* - 252(/i - lf{u - - 657) 



287(m - 



30375 

+ 3(m 1)^(50647^* - 43939) - 2{u - l)-^(29858/i - 23191)] 

- 48(5(At(2z/ + 21) - lv{ii - 7)) [6/^^ + ^^(Jy - 1) + 7!/(4!/ - 1)] ^ = , (3.11) 

40320a4!^-^ [6/i-'' + 'iii^iSv - 7) + lA^iv{y - I) + {v - l)y^] 

- '^^^^'^^^^ [l7739(/i - 1)^ - 27(At - 1)^(15259;^ - 12631) + 13334(i/ - lf{v'^ - 9zv + 5) 

- 9(/7, - l)-''(94639i^2 - 102980;/+ 20167) + 3(m - - l)(362107iy2 - 518030^+ 162475) 

- 2(^ - - 1)2(258865^^ - 468650^ + 179599)] 

- 165(6/^2 + 47m + 7(At + 25)z/ - 281/^ - 168) [6/1^ + ^iiilv - 1) + ly{Ay - 1)] ^ = , (3.12) 

17280a4Mi^'' [5a*^(3/x - 7) + 2iiu{23ii - 21) + (38At - 7)u'^ + 6iy^] 

+ ^^3Q3~g^' [41391(M - 1)' - 9(M - 1)^(70015;. - 58189) 

- 3(m - 1)3(151793j/2 - 58702i/ - 57613) + (m - ^fi^ - l)(784721zy2 - 937318i/ + 172253) 

- 2{fj, - l)(zy - l)2(202339z/2 - 338869i/ + 106344) - 40002(z/ - lf{2u - 1)] 

+ 48(5(/x(2m + 3z/ - 25) - 49z/ + 56) [6/^^ + 3^(7,^ - 1) + 7^(4^ - 1)] ^ = . (3.13) 
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Two of them are independent as we have shown. Using the values for 0:4 and 7 in Eq. (2.26), we have solved 
these equations numerically and found the following ten solutions for ^ = 0: 

{11, ly) = MP1(1.58841, 0.319 25) , MP2(0.733 61, 0.082 88) , MP3(0.722 46, -0.166 74) , 

MP4(0.622 07, -0.40040) , MP5(0.100 32, -1.700 96) , MP6(0.022 04, 0.990 55) , 

MP7(-0.030 14, 0.620 90) , MP8(-0.335 30, 0.850 24) , MP9(-0.668 48, 0.634 27) , 
MP10(-0.938 01, 2.572 50) . (3.14) 

3.2 (7p = 0, CTg 7^ (or ap j^O,ag = 0) 

3.2.1 Generalized de Sitter solutions (e = 0) 

Here we have Ap = ij?, Aq = + aqe~'^'^*, X = and Y = v"^, where dq = aq/bQ. It is easy to see that there 
is no exact solution unless = 0, in which case we have constant Ap = X = /j,^, Aq = dq and Y = 0. Our 
basic Eqs. (2.17) and (2.19) now read 

ai [pii/ + qiaq] + m \p-ji^ + ^PsqilJ^^^q + Qpsqsli^^q + ^piq^lJi^dq^ + q^dq^] 

+ 5[p{p-7)l? +ql^q][pl^J? +qidqf =Q , (3.15) 

ai[(p+ l)oM^ + {q- 1)2^9] 

+ a4 [{p + 1)6/ + 4(p + 1)492/^5 + 6{p + 1)294/^^^2 + 4(p + l)o%M^^g + qs^t] 

7(g- l)(p+l)ojVi(g,p) ^ ^ ..nr2,~]3 
^ (£)- l)3(j_2)4 +(9-8)a,J[M +aq\ 

+ 6[ip+ l)o/x2 + {q- l){q - 8)aJ [{p + 1)om' + Qi^g] ' = . (3.16) 

We note that Eq. (2.18) gives the same equation as Eq. (2.19) for = Q and need not be taken into account. 
Setting (5 = 0, it turns out that there is no solution for ap = 0,aq ^ 0. For the case of fjp ^ and aq = 0, 
exchanging and v, q, we find that there is also no solution. Although there is no this kind of solution for 
(5 = 0, we find solutions for ^ ^ as discussed in § 4.2.1. 

3.2.2 Pow^er-law solutions (e — 1) 

Here we have Ap = ii^,Aq = v"^ + a^e^*^^"^)*, X = /z(/i — 1) and Y = v{v — 1). We have only asymptotic 
solutions in most cases. 

{l)v>l: 

For t — > 00, the EH term dominates and Aq v^. The solutions are the same as C7p = (Tg = case in 
section 3.1.2. However, there is no solution with v > \. 

For t — > —00, Aq aqe'^^^~''^* and there is no solution. 
(2) u<l: 

For t — » 00, Aq (TqC^'^^^''^* and there is no solution. 

For t —00, Aq and the solutions are the same as (Tp = (T, = case in section 3.1.2. Choosing 

solutions with v < \ from Eq. (3.14), we get nine solutions MPl - MP9. 

We have Ap — , Aq = \ ^ aq,X = /i(/i - 1) and y = 0. 
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For t ^ CO, the EH term dommates and the solution is 

(At, u, dq) = ME12(0, 1, -1) . (3.17) 

Actually this is an exact solution. 

For t — > — oo, fourth-order terms dominate. Our basic independent Eqs. (2.17) and (2.19) give 

cti [p7^J^^ + Speg/i^ + 4p5(7iAi^(6 + Aq) + 8^492^^(4 + 'iAq) + 2p^qs.^i.'^{S + 24A, + SAg^) 
+ 8p2g4A*^^g(4 + iAq) + 4pi(?5Ai^Aq2(6 + Aq) + 8pqenAq^ + 97^/] 

+ iD-%{D-2r ~ ^^"^ ^ ^^^^ ^ '^^ ~ ^^^^ - 8(P + 9 - 6) + (P + 1) A] [m(m-2) + AJ' 

+ <5[p(p - 7)/x' + 2p(g - 21)m - 48? + q,Aq\ [p{p + l)^^ + 2p{q - 1)^ + giAj^ = , (3.18) 



"4 [{P + 1)6M* + 8P5(M ~2p + 6)m^ + 4p4(9 - 1)m'(6(p9 - 4p - g + 12) + (p + l)(g - 2)^,) 
+ 8p3(g - l)2M'(4p(g - 6) - 8(g - 9) + 3{pq -~Ap + 4)A,) 

+ 2p2(g - 1)3m'(8(p - 3)(9 - 8) + 2A{pq ~ 6p - q + 10) Aq + 3(p + l)(g - 4)A/) 
+ 8pi{q - l)4M'^g(4(p - 2)(g - 8) + 3(m - 6p + 2)^,) 

+ 4p(g - 1)5m'A'(6(p - l)(g - 8) + (p + l)(g - 6)A,) + 8p(q - 8)(g - 1)5^^^/ + {q - 1)3^/] 

+ g-^)3p-2^^ [^^^ - + ' + - + - + 
-8{q-8)ip + q-6) + ip+ l){q - 8)Ag)] 

+ S[(p+ 1)om' + 2p{q - 29)^i + {q - 8){{q - l)Aq - 48)] [{p + 1)om' + Mq ' 1)m + 9i A] ' = • (3-19) 

For p = 3 and q = 7,we find the solutions MEIO in Eq. (3.17) and 

(/i,i/,CTq) = MP11(6.679 58, 1 -2.106 62) , MP12(6.085 83, 1, -1.35793) , MP13(2, 1, -1) , 

MP14(0.481 80, 1, -1.238 02) , MP15(0.072 89, 1, -2.208 58) . (3.20) 

For the case of cTp ^ and aq = 0, exchanging and we obtain one exact solution ME13 and 
thirteen asymptotic solutions MP2 - MPIO in Eq. (3.14) and the following solutions MP16 - MP19: 

(/i,!/,ap) = ME13(1, 0, -1) , (3.21) 
(/i,!/,CTp) = MP16(1, 0.26818, 1.43812) , MP17(1, -9.17779, 4.633 79) , 

MP18(1, 0.93176, -3.96511) , MP19(1, -0.122 50, -1.414 71) , (3.22) 

where dp = ap/oQ. 

3.3 CTpCTg ^ 

3.3.1 Generalized de Sitter solutions (e = 0) 

li jjL = V = Q, our basic Eqs. (2.17) and (2.18) reduce to 

Oil [Pi^p + li^q] + ai[p7dp'^ + Ap^qiCFp^dq + &p?,q?,d/dq^ + Apiq^apdq^ + q7dq'^] 
^ {D - 1)3(£) - 2)4 [^'^^^ + 1)0^1 (P> <l)^P^ + 4pigiiVi(p, q)dp^dq 
+ 2piqiN2{p, q)crp^dq^ + ipiqiNi{q,p)dpdq^ + {p + 'i-)oqiNi{q,p)dg^] 

+ S[pidp + qiaqY = , (3.23) 
ai [{P - l)2CTp + qidg] 
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(3.24) 



+ 2p,(q - 4)(rj - 1)7V2(P, g)5p'<^,' + 4pi((z - 6)(g - l)iVi(g,:p)apa/ + (p + l)o(g - 8)(g - l)iVi(g,p)a/] 



For p = ?i. q — 7 and 5 = 0, wc find that there is no solution. 

If /U ^ and = 0, it is clear that there is no exact solution. For asymptotic solutions, we can search for 
them by setting Ap = ji^, Aq = dq, X = ji^ and F = 0. This is the same condition in § 3.2.1, and we have no 
solution. The case with = and ^ is obtained by exchanging p, fi and q, v. We find again that there 
is no solution. 

For fj,!^ ^ 0, if our ansatz for solution is imposed, it is easy to see that there is no solution if fi and v are 
of the opposite signs. If they are of the same sign, either t +oo ov t ^ —oo gives Ap ^ jj?, Aq ^ and 
there may be solutions. For 5 = 0, however, we see that Eq. (3.4) gives no solution of the same sign. 

3.3.2 Power-law solutions (e = 1) 

In this case, we first consider the cases when both [i and v are not equal to 1. 

(1) ii>\ and 1/ > 1: 

For i — !■ oo, the EH term dominates and we obtain the asymptotic solutions in § 3.1.2. Again no solutions 
satisfy the condition of /x > 1 and v > \ (see Eq. (3.8)) and hence there is no asymptotic solution of our 
form. 

As t — > — oo, the fourth-order terms become dominant and we find no consistent solution from the 
forth-order terms. 

(2) < 1 and v <\: 

As t ^ oo with EH dominance, wc again find no consistent solution. As t ^ — oo with fourth-order-term 
dominance, we obtain eight asymptotic solutions MP2 - MP9 from Eq. (3.14) in § 3.1.2: 

(3) \x>\ and v <\: 

As f ^ 00, ^ and Aq — ^ aq^'^^~"'>^ . This is similar to the case (2) in § 3.2.2 and there is no 
solution of our form. 

As t ^ -oo, Ap CTpe^^^"'*^* and Aq . We find no solution. 

(4) n< \ and v > \\ 

Hero wc reach the same result by exchanging p, /x and q, v. No asymptotic solution of our form is obtained 
for both t ±oo. 

Next, we discuss the cases in which one of or i/ is equal to 1 and the other is not: 



+ ^[(PiCTp + (g - 8)(g - 1)0-5) (piCTp + 9i^g)^] = 



(3.25) 



where q) is defined by 

A^2(P, 9) = P^(99 - 16) +p^(3g2 - 1) + ^^(M?^ _ 2\q + 26) + Zv^{q - 3) - 12/g3 

- 13p2g2 + + _ 3)^2 ^ _ + 26)g3 -h (3^^ - l)g4 + (g^ _ i6)g5 . 



(3.26) 
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(5) > 1 and ly — 1: 

As f — > 00 with EH dominance, Ap y?, and we recover the case of cTp = 0, ^ 0. However, there 
is no solution with /i > 1. Wc do not liavo any asymptotic solution of our form. As t ^ — oo with 
fourth-order-term dominance, Aj, o-pC^'^"'^'*. We again do not have any asymptotic solution of our 
form. 

(6) ii<l and v = 1: 

As f — > 00 with EH dominance, Ap diverges as tTpe^^-^"''^*. There is no asymptotic solution of our form. 
As t ^ — oo, we again recover the case of ct^ = 0, (Tq 7^ with the fourth-order-term dominance and 
solutions in (3.17) and (3.20). (Note that (3.17) was an exact solution for ap = 0.) Choosing those with 
< 1, we get asymptotic power-law solutions 

(/i,i^,a,) =MP23(0, 1, -1) . (3.27) 

form Eq. (3.17), and MP14 and MP15 from Eq. (3.20). 

(7) /X = 1 and z/ > 1: 

The analysis is almost the same as the case (5). There is no asymptotic solution. 

(8) /i = 1 and v <l: 

The analysis is almost the same as the case (6), and we find the asymptotic solutions 

(/X, V, ag) = MP24(1, 0, -1) . (3.28) 

form Eq. (3.21), and MP16 - MP19 for t — > —00, which are the same as the case of ap ^ 0,aq = given 
in Eq. (3.22). 

Finally, we consider the remaining case. 

(9) /X = 1 and v = 1: 

Here we have constant Ap = 1 + ap and Aq = 1 + aq. As t ^ +00, the EH term is dominant, and we 



have 



PiAp + qiAg + 2pq = 0, 

{p - l)2Ap + qiAq + 2{p - l)q = 0, 

PiAp + {q- 1)2^, + 2p{q - 1) = 0. (3.29) 



The solution is given by 



p~V q-V 



Ap = 2^, Aq^--^. (3.30) 



This is the solution found in Ref. [11] which exhibits eternal accelerating expansion when higher order 
effects are taken into account. For p = 3,q = 7, we get the following future asymptotic solution: 

{li,iy,ap,aq)=MF8(l, 1,-1,-1) ■ (3-31) 



2' 2^ 

For t ^ —CO with fourth-order- term dominance, we get two independent equations for dp and aq from 
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Eqs. (2.17) and (2.18). 

ai [p7{l + dpf + 4(1 + apf{2qp(i + p^qi{l + 5-,)) + 24(1 + apfipsQi + ^492(1 + a,)) 
+ 32^492(1 + CTp) + 32p294(l + ^q) + 24(pig5 +p294(l + ffp))(l + ^qf 
+ 4{2pqe + Piq^il + ap))(l + agf + q^il + agf + 6p3(Z3(l + ^pf{^ + ^gf 
+ 48p393(l + ^p)(l + + I6P393] 
+ p _ -^^^3^^J) _ [P^il + 1)0^1 (P, q)(^p^ + 4:piqiNi{p, q)dp-^CFq + 2piqiN2{p, q)cFp^dq^ 

+4piqiNi{q,p)apaq^ + {p + l)aqiNi{q,p)aq^] 
+ S[{{p + q-49){p + q)+Pidp + qiag){{p + q-l){p + q)+piap + qiagf] = , (3.32) 

"4 [ip - 1)8(1 + ^pf + 4(1 + ^pf{2q(p - 1)7 + (p- 1)6(71 (1 + ag)) 

+ 32{p - 1)5^2(1 + ap) + 32{p - 1)3^4(1 + ag) + 24{{p - 1)2^5 + {p - 1)3^4(1 + ap)){l + dgf 
+ A{2{p - 1)96 + {p- 1)2^5(1 + <?p))(l + ^g)^ + 97(1 + ^g)* + 6(p - 1)493(1 + ^pf{^ + ^qf 
+ 48(p - 1)4^3(1 + ^p)(l + dg) + 16(p - 1)4(73] 

+ (£) - lf{D - 2)4 ~ - 8)(9 + l)o^p + 4(p - 6)giag - 8(p + g - 7)^} 

- qiNi{q,p)ag^{A{p - l)2ap + {jp + l)oag + 8(p + g - 7)^} 

+ 4(p - l)q^{p + q- 7fapag{N3{p, q)ap - N3{q,p)ag) - 2{p - l){p - 4)7V2(p, q)^p'^^q''] 

+ S[{ip + q- 49)(p + g -8) + {p- 8){p - l)ap + qiag){{p + q- l){p + q) +prap + q^dgf] = , 

(3.33) 

where we define N^{p, q) as 

N^iv, q) = Wi( V + 35p - 18) - 3pi(6p2 -p-&)q- (32p3 - 4^ + 33j3 - 18)5^ 

+ (16p2 + I5p - Al)q^ + (21p - 41)q'' - 2q^ . (3.34) 

For p = 3, q = 7 and S — 0, we find four solutions 

{li,u,ap,ag) =MP25(1, 1, 2.37016,-1.262 93) , MP26(1, 1, -1.412 83,-120.598 08) , 

MP27(1, 1, -4.376 02,-0.120 21) , MP28(1, 1, -5.762 67,-2.588 03) . (3.35) 

We summarize our results obtained here in Appendix D. The exact solutions for (5 = are listed in Table 5, 
future asymptotic solutions in Table 6 and past asymptotic solutions in Table 7. The numbering of solutions 
are given for those with ^ ^ (see next section). For example, in Table 5, we do not find solutions ME2± - 
ME11±. 

4 Solutions in M theory for 6^0 

In this section, we search for exact generalized de Sitter solutions for —10 < (5 < 10 whereas we give a few 
examples for future and past asymptotic solutions for particular 5 to avoid vexatious complications. As typical 
examples, we mainly focus on the case oi S ~ —0.001 and S = —0.1 as to future and past asymptotic solutions 
in § 4.1.2, 4.2.2 and 4.3.1 - 4.3.2, but we also study how the results change depending on 6. 

4.1 ap — aq — 

4.1.1 Generalized de Sitter Solutions (e = 0) 

In this case, we have the same Eqs. (3.1) and (3.2) with a free parameter 6^0. For the given value of 6, we may 
solve these equations numerically. In Figs. 1 and 2, we depict numerical solutions MEi+{6, fJ,i,yi) (i = 1, • • • , 5) 
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with fii > when 5 is varied. We note that there are always time- reversed sohitions MEi-{6, fi'^, i/'^) {i = 
1, • • • ,5) obtained by (/i^jt'j') = {—jjii^—Vi) which are not shown exphcitly. We find five solutions ME1+ - 
ME5+ for 5 < 0, while just one ME1+ for J > 0. For the case of 5 = 0, we have a solution ME1+ which is 
consistent with the result of § 3.1.1. In Table 1, we summarize their properties. 
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Figure 1: Five Generalized de Sitter solutions with ap = Figure 2: One Generalized de Sitter solution with ap = 
CTg = with respect to 5 < 0. a, = near the origin of S. 

Each pair of points {5,fj,i,Vi) with the same value of 5 gives one solution. Number of solutions changes with 
the value of 5. There is another set of time-reversed solutions MEi_ with /ii < 0. No solution exists for 10~^ < 5. 



Table 1: Generalized de Sitter Solutions MEi+ (i = 1, • ■ ■ ,5) with fii > for various values of S. Five eigenmodes for 
linear perturbations are also shown. (ms,nu) means that there are m stable modes and n unstable modes. The solution 
has many stable modes if its 10- volume expansion rate 3/x -I- is positive. 



Solution Property Range Stability 3/ii + Ivi 



ME1+ 


lyi < < fii 


-0.000 04311 <S <0 


(0s,5u) 








5 = 


(ls,2u) 








0<5 < 0.000 059 88 


(ls,4u) 




ME2+ 


1^2 <0 < tJ-2 


-0.045 20 <S < -0.002 649 


(4s,lu) 


+ 


ME3+ 


< 113 = 1^3 


S < -0.000 4732 


(3s,0u) 


+ 


ME4+ 


V4 < < IJ.4 


-0.2073 <5 < -0.004 852 


(ls,4u) 




ME5+ 


i'S < < ^5 


-0.2073 <S < -0.2056 


(2s,3u) 





We find that there are two solutions around the value S —0.001 in Fig. 1. Especially, we have the 
following solutions for S = —0.001: 

{S,n,iy) = ME2+(-0.001, 1.43787, -0.0670662) , ME3+(-0.001, 0.402934, 0.402934) . (4.1) 

It is interesting to see how the solutions change for other value of 6. For instance, we have the following two 
solutions for S = —0.1: 

{S, fi, v) = ME3+(-0.1, 0.168 203, 0.168 203) , ME4+(-0.1, 0.742 918, -0.453 997) . (4.2) 
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where we use the same names for the sohitions connected when 6 is changed. The value of /x, and Vi of the 
solution with the same sign change with the value of S as in Eqs. (4.1) and (4.2). We will also study a linear 
perturbation around these solutions in § 5.1. 

4.1.2 Power-law solutions (e = 1) 

There is no exact solutions, but we have asymptotic solutions for various values of S. We will give them 

explicitly for S = -0.001 and 5 = -0.1. 

As t — » 00 with EH dominance, we get the same solutions MF6 and MF7 given by Eq. (3.10) in § 3.1.2 
As t ^ — oo, the forth order terms dominate. We have the same Eqs. (3.11) - (3.13) in § 3.1.2 with 6^0. 

For (5 = —0.001, we find the following twelve solutions 

((5,/i,j/) = MP6(-0.001, 121.218, -5.48783) , 
MP8(-0.001, 26.6578, -37.1453) , 
MP10(-0.001, 0.737 553, -0.086 3059) , 
MP12(-0.001, 0.190 928, 0.139 524) , 
MP14(-0.001, 0.120 104, 0.169 926) , 
MP16(-0.001, -1.161 61, 1.49783) , 

while the following eight for S = —0.1 

{5,iJ,,u) = MP6(-0.1, 14.0692, 14.0692) , 

MP8(-0.1, 0.229 099, 0.151 155) , 
MP10(-0.1, 0.123 489, 0.196 532) , 
MP12(-0.1, -60.445, 36.7782) , 

4.2 (Tp = 0, CTq 7^ (or (7p ^ 0, (Tg = 0) 

4.2.1 Generalized de Sitter solutions (e = 0) 

Here we have Ap = fj?', Ag = i/"^ + o-^e"^''*, X = /z^ and Y = v'^. It is easy to see that there is no exact solution 
unless = 0, in which case wc have constant Ap = X = , Aq ^ and F = 0. Our basic equations reduce to 
Eqs. (3.15) and (3.16). In Figs. 3 and 4, we depict numerical solutions MEi_|_(J, /Zj, t-j = 0,CTg) (i = 6,7) with 
/Uj > as a function of b. There are also time-reversed solutions MEi_((5, /U-, z^i = 0, {i = 6, 7) obtained by 
/i^ = —/;,,; which arc not shown. We find a solution ME6-|- for S < 0, while ME7+ for 5 < 0. In the vicinity 
of (5 = 0, we find no solution as discussed in § 3.2.1. In the anterior part of Table 2, we summarize their 
properties. 

Table 2: Generalized de Sitter Solutions MEi+{5, /lijUi) (i — 6, . . . , 11) with j-Li,i^i > for various values of 5. Six 
eigenmodes for linear perturbations are also shown. (ms,nu) means that there are in stable modes and n unstable 
modes. 



Solution 




Property 


Range 


Stability 


3/ii + 7i>i 


ME6+ 


fe = 


= 0, < /J6, Orq(6) 


5 < -0.000 5589 


(5s,lu) 


+ 


ME7+ 


vy = 


0, 5-g(7) <0 < He 


0.002 999 < S 


(4s,2u) 


+ 


ME8+ 


fis -- 


= 0, < 5-p(8), 


-0.003 163 < (5 < -0.000 5650 


(4s,2u) 


+ 


ME9+ 


M9 -- 


= 0, < ap(g), i/g 


5 < -0.000 5657 


(5s,lu) 


+ 


ME 10+ 


Mio = 


0, 5-p(io) < < uio 


S < -0.000 043 49 


(5s,lu) 


+ 


ME11+ 


Atii = 


0, 5-p(ii) < < 1^11 


-0.085 22 <S < -0.003164 


(4s,2u) 


+ 



Especially, we find the following solution for S = —0.001: 

((5, At, u, dg) = ME6+(-0.001, 0.775 074, 0, 0.278 981) . (4.5) 



MP7(-0.001, 27.0789, 27.0789) , 
MP9(-0.001, 2.61038, -0.118736) , 
MP11(-0.001, 0.726 753, -0.151 41) , 
MP13(-0.001, 0.154 834, 0.154834) , 
MP15(-0.001, -0.757 551, 0.625 032) , 
MP17(-0.001, -2.407 56, 0.598 625) , (4.3) 

MP7(-0.1, 8.240 22, -2.53152) , 
MP9(-0.1, 0.174 972, 0.174972) , 

MP11(-0.1, -5.57678, 3.39005) , 
MP13(-0.1, -225.859, 69.5863) . (4.4) 
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Figure 3: One generalized de Sitter solution with (Tp 
0, o-q for J < 0. 
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Figure 4: One Generalized de Sitter solution with cTp 
0, for < 5. 



while the following for 6 = —0.1 

(5, [1, V, dq) = ME6+(-0.1, 0.307 198, 0, 0.047 1560) , 



(4.6) 



For the case of o-p 7^ and aq ~ 0, exchanging /x, p and z/, g, we find four solutions ME8+ - ME11+. In 
the vicinity of 5 = 0, we find again no solution as discussed in § 3.2.1. In Fig. 5, we depict these four numerical 
solutions with Vi > with respect to S. There are also time-reversed solutions MEi_((5, fii = 0, v^, dp) (i = 6, 7) 
obtained by i>l — —vi which are not shown. In the posterior part of Table 2, we summarize their properties. 

Especially, for 5 = —0.001 and 5 = —0.1, we have the following solutions, respectively: 

((5, /X, ly, CTp) = ME8+(-0.001, 0, 0.491 829, 0.911 910) , ME9+(-0.001, 0, 0.401 567, 1.970 26) , 

ME10+(-0.001, 0, 0.408 387, -0.672 260) . (4.7) 
ME9+(-0.1, 0, 0.201057, 0.141557) , ME10+(-0.1, 0, 0.314065, -0.757016) . (4.8) 

4.2.2 Power-law solutions (e = 1) 

Here we have Ap = 11^, Aq — + (TgC^^^""^'*, X = /x(/x — 1) and Y = v{v — 1). We have only asymptotic 
solutions in most cases. 

(1) i^>l: 

In this case, we have same result in § 3.2.2 (1), and there is no asymptotic solution of our form. 

(2) V <l: 

For t 00, Aq — > a-qC^^^^'^^* and there is no solution. 

For t —00, Aq and the solutions are the same as o-p = cr, = case in § 4.1.2. From Eqs. (4.3) 

and (4.4), we get ten past asymptotic solutions MP6, MP8 - MP15 and MP17 for 6 = -0.001, while the 
four solutions MP7 - MPIO for 8 = -0.1. 



(3) = 1 : 
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Figure 5: Four Generalized de Sitter solutions with Op ^ Q, aq = G with respect to (5 > 0. There is no solution for 
<5 > 



We have Ap ^ fM^^Aq ^ 1 + ct,, X = - 1) and F = 0. 
For t — > (X), the EH term dominates and the solution is 

=ME12(0, 1, -1) 



(4.9) 



Actually this is an exact solution for all S. 

For t — oo, fourth-order terms dominate. Our basic Eqs. (2.17) and (2.19) give the same Eqs. (3.18) 
and (3.19) in § 3.2.2 with S 0. We get the following two solutions for 6 = —0.001 

id,n,iy,ag) = MP18(-0.001, 32.4979, 1, 482.327) , MP19(-0.001, -47.5769, 1, -172.534) , (4.10) 

while the following for 6 = —0.1 

{S,fi,iy,dg) = MP14(-0.1, 14.0811, 1, 92.0567) , MP15(-0.1, -31.2288, 1, -201.029) , (4.11) 

For the case of dp 7^ and aq — 0, exchanging and v^q, we obtain one exact solution ME13 and 
eleven past asymptotic solutions MPIO - MP17 in Eq. (4.3) and the following solutions MP20 - MP 22 for 
S = —0.001, while one exact solution ME13 and nine past asymptotic solutions MPS - MP13 in Eq. (4.4) and 
the following solutions MP16 - MP18 for S = -0.1: 



(m, i^, (Tp, ao) = ME13(1, 0,-1,1) . 
{S,fi,iy,ap) = MP20(-0.001, 1, 31.7651, 3563.36) , 

MP22(-0.001, 1, -0.84176, -1.25634) , 
MP16(-0.1, 1, 14.0787, 618.521) , 
MP18(-0.1, 1, -279.628, -1.00255 x 10*^) 



(4.12) 

MP21(-0.001, 1, 24.9773, 8412.47) , 

(4.13) 

MP17(-0.1, 1, -0.634 680, -1.447 37) , 

(4.14) 
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4.3 CTpCTg 7^ 

4.3.1 Generalized de Sitter solutions (e = 0) 

E H = v = 0, our basic equations reduce to Eqs. (3.23) - (3.25). We find again tliat tliere is no solution for 
^ ^ 0. If either = or i/ = and the other is nonzero, it is clear that there is no exact solution. For 
asymptotic solutions, we can search them by setting Ap = /j,'^ , Aq = aq, X ~ jj,'^ and Y = for the latter case. 
This case is actually the same as § 4.2.1, and thus Eqs. (4.5) and (4.6) give asymptotic solutions with = 0. 
For (5 = —0.001, we have the following future and past asymptotic solutions: 

{S, fj,, V, dq) = (-0.001, ±0.775 074, 0, 0.278 981) : MF2(MP2) for i ^ ±oo. (4.15) 

For 5 = —0.1, we have 

[8, n, V, dq) = (-0.1, ±0.307 198, 0, 0.047 1560) : MF2(MP2) for i ^ ±oo. (4.16) 

The first case is obtained by exchanging p, fi and q, u. From Eqs. (4.7) and (4.8), the solutions are 

{S,fi,u,ap) = (-0.001, 0, ±0.491829, 0.911910) : MF3(MP3) , 
(-0.001, 0, ±0.401567, 1.970 26) : MF4(MP4) , 

(-0.001, 0, ±0.408 387, -0.672 260) : MF5(MP5) for t ^ ±oo , (4.17) 

for 6 = -0.001 and 

{6,ii,v,dp) = (-0.1, 0, ±0.201057, 0.141557) : MF3(MP3) , 

(-0.1, 0, ±0.314065, -0.757016) : MF4(MP4) for i ^ ±oo , (4.18) 

for 6 = -0.1. Note that Eqs. (4.17) - (4.18) were exact solutions for o-p = or Uq = in § 4.2.1. 

For yu, 1/ ^ 0, if our ansatz for solution is imposed, it is easy to see that there is no asymptotic solution if 
/U and u are of the opposite signs. If they are of the same sign, either t +oo or t ^ —oo gives Ap ii^, 
Aq and there may be solutions which eventuate in these of § 4.1.1 with the same sign. This implies that 
the inflationary solutions with positive eigenvalues are obtained for asymptotic infinite future, so that they 
are not interesting from the cosmological point of view. However, it may be useful to check if there are any 
solutions of this type. In fact, Eq. (4.1) gives set of asymptotic solutions for b = —0.001 

(5, A«, i^) = (-0.001, ±0.402 934, ±0.402 934) : MFl(MPl) for t ^ ±oo . (4.19) 
meanwhile, Eq. (4.2) gives solutions for b = —0.1 as 

{S, v) = (-0.1, ±0.168 203, ±0.168 203) : MFl(MPl) for t ^ ±oo , (4.20) 
Note that Eqs. (4.19) and (4.20) were exact solutions for (jp,aq = in § 4.1.1 

4.3.2 Pow^er-law solutions (e — 1) 

In this case, we first consider the cases when both and u are not equal to 1. 

(1) /z > 1 and v>l: 

In this case, we have same result (1) in § 3.3.2, and there is no asymptotic solution of our form. 

(2) /X < 1 and v <l: 

As t ^ 00 with EH dominance, we again find no consistent solution. As t — > — oo with fourth-order- term 
dominance, and we have same result in § 4.1.2 with /x < 1 and v < 1. We obtain the seven asymptotic 
solutions MPIO - MP15 and MP17 from Eq. (4.3) for 5 = -0.001, while three asymptotic solutions MP8 
- MPIO from Eq. (4.4) for 6 = -0.1. 
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(3) iJ> 1 and i/ < 1: 

The analysis is almost the same as the case (3) in § 3.3.2, and there is no asymptotic solution of our 
form for both t ±oo. 

(4) /i < 1 and v > 1: 

The analysis is almost the same as the case (4) in § 3.3.2, and there is no asymptotic solution of our 
form for both t ±00. 

Next, we discuss the cases in which one of or 1/ is equal to 1 and the other is not: 

(5) > 1 and v = 1: 

The analysis is almost the same as the case (5) in § 3.3.2, and there is no asymptotic solution of our 

form for both t ±00. 

(6) fi < 1 and ly = 1: 

As t ^ 00 with EH dominance, Ap diverges as CTpe^*^^"'')*. There is no asymptotic solution of our 
form. As f ^ —00, we again recover the case of cTp = 0, fjq ^0 with the fourth-order- term dominance 
and solutions in Eqs. (4.9), (4.10) and (4.11). (Note that Eq. (4.9) was an exact solution for ap = 0.) 
Choosing those with /x < 1, we get asymptotic power-law solutions 

(/x,z^,aJ=MP23(0, 1, -1), (4.21) 

for all 5, and MP19 in Eq. (4.10) for S = -0.001 and MP15 in Eq. (4.11) for 5 = -0.1. 

(7) = 1 and v > 1: 

The analysis is almost the same as the case (5). There is no asymptotic solution. 

(8) = 1 and v <1: 

The analysis is almost the same as the case (6), then wc find the asymptotic solutions as t ^ —00, which 
are the same as the case of ap ^ 0,i7q = given in Eqs. (4.12), (4.13) and (4.14). We have asymptotic 
power-law solutions 

(Ai,i.,ap) = MP24(l, 0, -1) , (4.22) 
for all 6, and MP22 in Eq. (4.13) for 5 = -0.001 and MP17 and MP18 in Eq. (4.14) for S = -0.1. 
Finally, we consider the remaining case. 

(9) = 1 and v = 1: 

Here we have constant Ap = 1 + dp and Aq = 1 + dq. For t — > -|-oo, we have the same solution MF8 in 
§ 3.3.2. 

For t —00 with fourth-order-term dominance, we get same two independent Eqs. (3.32) and (3.33) in 
§ 3.3.2. We find solutions 

{S, dp, dq) = MP25(-0.001, 116.501, 9.455) , MP26(-0.001, 33.9247, 10.4469), 

MP27(-0.001, -0.886 667, -2.590 21) , MP28(-0.001, -3.020 49, -1.423 539) , (4.23) 

MP25(-0.1, 19.5872, 6.525 89) , MP26(-0.1, -219.014, 20.5567), 

MP27(-0.1, -0.448 085, -2.186 91) , MP28(-0.1, -3.443 20, -1.582 47) . (4.24) 

We summarize our results obtained in this section in Appendix D. Exact solutions for 5 = —0.001 are 
listed in Table 8, future asymptotic solutions for in Table 9 and past asymptotic solutions for in Table 10. For 
(5 = —0.1, we summarize only exact solutions in Table 11. 
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5 Stability Analysis of Generalized de Sitter Solutions 



Since the exact generalized de Sitter solutions MEl-t - MEll-t obtained in § 3.1.1, § 4.1.1 and § 4.2.1 correspond 
to fixed points in our dynamical system, we have to analyze their stabilities in order to see which solutions 
are more generic and also to find interesting cosmological solutions. Wc have performed a linear perturbation 
analysis around those fixed points for the solutions. In the following subsections, we classify the result by the 
signature of Gp and <jq. 

5.1 Up — (Jq — 

In this case, we have exact solutions ME1± - ME5± obtained in § 3.1.1 and § 4.1.1. Setting 

=Hi + Aie" * , -J- = Pi + Bie" * , (5.1) 
where \Bi\ <C 1, we write down the linear perturbation equations: 

GA{^^^,V^,^'''^)A,+GB{^l^,V^,^'''^)B,=0 , (5.2) 

obtained from Eqs. (2.17) and (2.18). Quantities Fa, Fb, Ga and Gb are functions of /Uj, Vi and cr^*-' given by 
Fa = 6ai [2^ + 7i^] + 120960a4i/^ [21/^^ + 14^i/ + i/^] 

- 7 [23652/ - 568368^*^1/ + 694341/^^^ - 301231/ii/^ + 29858i/* 
30375 '- 

- (16647/i^ - 8892/zi/ + 3313i/2)(3/u + 6i/ + a)a] 

- 1925(6At2 + 2lnv + 2Sv^f [24// - 21/^^;/ - ll^jj.v'^ - 49u^ 

-3(9/|2 + 49Ait/ + 42i/2)o--3(3/i + 7i/)cr2] , (5.3) 
Fb = 42ai [fi + 2u] + 846720a4Ati/^ (a* + v) [5^J, + u] 

+ ^^^30375^''' 7[At(91332/x^ - hAMll^i^v + 545397mz/2 - 209006z/3) 

- (16647/i^ - 8892/iz/ + 3313;/^) (2/* + 7v + u)u] 

- 1344(5(6^^ + 2liiv + 2Sv'^f [^(3^^ - 9^;^ - 491/^) - (6At^ + 45;uz/ + A'dv^)a - (3;U + 7v)a'^] , (5.4) 
= 4ai [3/i + 7t/ + cr] + 8G64Ga4i/^(3/x + v) [3/x + 7i/ + ct] 

- 567(/^ - '^^ '^^ [23652/ - 550629^1/ + 707712/1/^ - 320185/iz/^ + 39838z/^ 

- (16647/ - 8892/iz/ + 3313z/^)(3m + 1v + u)a\ 

- 64(5(6/ + 2l^iv + 28i/^)^(3/i + 7i/ + a) [24/ + 21nv - 56i/^ - 9(3/x + 7;/ + C7)c7] , (5.5) 
Gb = 14ai [2/i + 8i/ + ct] + 40320q;4Z/'^(15/ + 12iiv + u^) [2/i + 8i/ + cr] 

+ ^^^Q^^25''^ [232605/ + 27765/z/ - 3758859/i/^ + 4294423/z/^ - 1845390Atz/* + 53336i/^ 

- (8550/ + 1094319/z/ + 80985/i/^ - 156303^2/^ + 178861z/^)cr 

- (16647/ - 8892/iz/ + 3313z/^)(5m + 15z/ + (t)c7^] 

- 448(5(6/ + 21/iz/ + 28z/^)^ [15/ + 45/z/ - 56/x!/^ - 22Av^ - (15/ + lll^i/ + 196j/^)cr 

- 3(5m + 15z/ + cr)cr^] , (5.6) 

for p = 3, q = 7, where we omit the subscript of in, and (T^*^. The condition that Eq. (5.2) has nontrivial 
solutions for Aj and Bi is 

i^A(Mi,z/i,aW)GB(/Xi,i/i,a«)-i^B(Mi,J^i,«^^'^)GA(Mi,i^i,«^^'^) = , (5.7) 
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which yields five modes (cr ~ aa\ a — 1, 2, • • • ,5) for each solutions i = 1, • • • , 5 with fixed S 0. This is 
because the basic equations for iii and U2 are two third-order differential equations plus one constraint which 
is second order. Eqs. (5.2) are derived from Eqs. (2.17) and (2.18), but we have checked that the results for 
(Ta ^ remain the same if we use any two combinations of Eqs. (2.17) and (2.19). For instance, we have five 
modes 

ME2± : cr(*) = (T5.364 98, T3.84414, ^3.773 01, tO.071 1343, ±1.520 84) , 

ME3± : o-W = (IF4.029 34, ^3.939 81, tO.089 5274) (5.8) 

for S = -0.001 and 

ME3± : cr(^) = (Tl.682 03, Tl-609 88, tO.072 1502) , 

ME4± : o-^') = (t3.82 81, ±0.072 7295, ±0.876 493, ±0.949 222, ±4.777 32) , (5.9) 
for S = —0.1. The class of solutions ME3± has only three modes because these solutions are special case with 

/U = !/. 

For S = 0, we have only three modes a^^\ (a = 1,2,3) because of the conformal invariance of the Weyl 
tensor. Specifically, we have the following three modes for ME1± 

ME1± : o-W = (IF3.871 09, ±6.242 68, ±10.1138) . (5.10) 

The numbers of stable and unstable modes for various values of S are summarized in Table 1. The number 
of unstable modes is important to discuss generality of inflation. For example, for the solution ME4_|_((5, H4,V4), 
there are one stable and four unstable modes. This implies that this solution may not be generic because there 
are many unstable modes. The probability to approach such generalized de Sitter solution will be very low. 
On the other hand, the solutionME2+ has four stable modes as well as one unstable mode. Hence, except 
for one direction in the phase space, this solution is stable. There may be a finite probability that a generic 
spacetime first approaches to this solution and eventually evolves into other solution. The solution ME3+ 
has three stable modes, which means that this solution is stable against linear perturbations. We find that 
preferable solutions, that is, the solutions with many stable modes are obtained when its 10-volume expansion 
rate (3/ij + 7iy'i) is positive. 

5.2 ap — 0,aq (or ap 0,aq — 0) 

We have exact solutions ME6± and ME7± obtained in § 4.2.1 for (jp = and aq ^ 0. We have also carried 
out the linear perturbation similar to that in § 5.1. Writing 

at 

= ln6o + Bie'''''*, (5.11) 

we derive the linear perturbation equations (5.2). We find 

Fa = 2piaiii + 8a4fi[p7lJ.^ + Spsqifi'^ag + Sp^qan'^aq'^ + piq^dg^] 

+ ^^^^^^^^^^[2m{q,p){{p-7)f,' + ip-l)aq) + {2N2{p,q) - {j, - m^{Q,P))iP^^ + M 

+ 8p5iJL{{p + l)oi/ + qidqf [(p -7){p+ l)oA<^ ± Piqiag + 6(pV + '^)<^] , (5-12) 
Fb = 2ai [pqiJ.a - q-^ag] - Sa^ [a-,(p59iM^ + ^PaQsl^'^^q + ^PiQsl^^^q^ + Qr^q^) 
- aipeqii^ + 3p4q2li^a-g + 3p2g4A*^a-,^ + pqej^ag^)] 

, 4pqai^i^ + agf ^2iVi(5,p){a,((5 - p)^' - (p + l)a,) + (g/x^ + {q- 6)/xa,)a} 



{D-lf{D-2Y 
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+ {2iV2(p, q)-{p- ?,)N:i{q,p)}nG{ii - a){py. + a)] 
+ 85((p + l)o/x^ + qic^qf [6(7^(7^ (pCT + pi/z) + pna{p(j> - 5)qii^ + (q - 6)qi^q} 

- - 5)a«^ + giCTj] , (5.13) 
Ga = 2pai [{p + 1)m + 0-] 

+ 8a4((p + l)At + ct) [pe/ + 3p4(5 - l)2M^CTg + ^P2{q - iW^^q^ + p{q - 1)6^9^] 
4p(g - l )j{{p+l)^i + a){n ^ + dqf 

-a{a+ p^)i2N2{p, q) - {p - i)Ni{q,p))\ 
+ ^p6{{p + l)n + (7)((p + 1)om2 + ^g^^p^ + ^) + (p + i)g^2 ^ _ _ g^~^j ^ 

Gs = 2ai [{q - l)u{pii + <j)-{q- 1)20-9] 

+ 8ai[a{piJ. + a){p5{q - l)lJ,^ + 3p3{q - 1)3/"^^^, + 5p2{q - 1)5/"^^^,^ + {q - l)?^,^) 

- aq{3{p + l)o(g - 1)6M^^,^ + (9 - l)8CTg^ + (9 - l)4M^((p + 1)3/"' + 3p2(p^ - P - 2)a,)}] 

- ^f^ll)Ip^2y [2iVi(,,p){(p + l)aqi{q - 2)^^' + {q- 8)a,) 
+ cr(p/i + (T)(2(2p + q- 2)fi^ + {p + 2q- 13)5-,)} 

- a(pn + a){2N2ip, q)a{pn + a) - N3{q,p){{l ~ p^)^i^ + {p - i)a{p^i + a))}] 

+ S5{{p+l)o|l^ + qlaqf[Qqa^{p^JL + af+a{p|l + a){{q-l){p+l)o^l^ + ql{q-V6)aq} 

-{p+ lUq - l)2lJ^aq -qi{q- 8){q - l)aq^] (5.15) 

The condition for the existence of nontrivial solutions (5.7) yields six modes (ct = aa \ a = 1, 2, • • • , 6) for each 
solution with fixed 6 because we have the new variable ui in addition to the variables in § 5.1. For instance, 
we have five modes 

ME6± : cjW = (IF3.305 74, ^2.177 00, Tl-162 61 T 1-516 84i, 

T 1.162 61 ± 1.516 84i, t0.148 218, ±0.980 516) , (5.16) 

for 6 = -0.001, while 

ME6± : a^^ = (t1.342 73, tO.770 998, t0.460 796 t 8.676 06i, 

T 0.460 796 ± 8.676 06i, tO.150 595, ±0.421 134) , (5.17) 

for 5 = -0.1. 

For cTp ^ and (Tq = 0, wc have four exact solutions ME9± - ME11± obtained in § 4.2.1. Exchanging fj,,p 

and ly, q, wc obtain the following six modes for each solution: 

ME8± : o-W = (T4.159 33, T3.348 59, Tl-721 40 T 2.744 09i, 

=F 1.721 40 ±2.744 09i, tO.094 2078, ±0.716 524) , 
ME9± : o-W = (=p3.ii5 29 T 1-707 36i, T3.115 29 ± 1.707 36i, T2.719 19, 

T 0.091 7844, ±0.304322 t 1-707 36i, ±0.304322 ± 1.707 36i) , 
ME10± : cr^'^ = (T3.30407, =f1-429 35 T 0.921 414i, =f1-429 35 ± 0.921 414i, 

T 1.429 35 T 2.731 49i, t1-429 35 ± 2.731 49i, ±0.445 363) for S = -0.001 , (5.18) 

ME9± : o-W = (iFl.733 58, Tl-31940, t0-703 698t 15.1267i, 

T 0.703 698 ± 15.1267i, t0.0879931, t0.326182) , 
ME10± : o-W = (T2.555 18, T2.408 32 T 1-076 03i , T2.408 32 ± 1.076 03i, 

± 0.209 864 T 1-076 03i, ±0.209 864 ± 1.076 03«, ±0.356 724) , for 5 = -0.1 . (5.19) 
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The numbers of stable and unstable modes for other values of S are summarized in Table 2. Solution ME6+ 
has five stable modes and one unstable mode, while ME7+ has four stable modes and two unstable modes. 



6 Solutions in Type-II Superstring 

We can discuss the case of type-II superstring in the same manner as M-theory if we keep the dilaton field 
constant and ignore the contributions from other fields. The low-energy effective actions for type-IIA and IIB 
superstrings with tree and one-loop corrections are [16] 



Sua = 
Sim = 

where an and 711 are given by 



1 



^10 



d^°xy/^ [R + a'^auEs + a'^juLw] , 
'g [R + a'^-fuLw] , 



au 



32-28 



711 



C(3) 
23 



2 2 

3-23 



(6.1) 
(6.2) 

(6.3) 



There is additional tree-level term Lw with coefficient ({3)/2^ in type-IIA superstring which vanishes under 
the de-compactification limit and does not exist in the M-theory action (2.1). The one-loop corrections in 
type-IIA superstring, i.e. Eg and the rest of Lw, are lifted to = 11 in the limit gs — )■ 00 and in agreement 
with M-theory corrections (2.5) and (2.6). We can again rescale an and 711 as Eq. (2.25), and have 



an 



•25 1^ ^n^g,^) 



711 



1 



(6.4) 



In the type-IIA superstring theory, we have a parameter < an < 1/(3 • 2^) corresponding to < g^ < +00, 
whereas an = in type IIB superstring. In the limit gs 00, the value of an is equivalent to a^ in M-theory. 
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Figure 6: Two generalized de Sitter solutions for type II 
superstring with cTp = cTq = with respect to < 04 < 
3~^2~^. a4 = corresponds to the type IIB superstring. 
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Figure 7: Two Generalized de Sitter solutions for type II 
superstring with o"p 7^ 0, ct, = with respect to < a4 < 
3~^2~^. Q4 = corresponds to the type IIB superstring. 



23 



We can search the exact and asymptotic sohitions in a same way in § 3, however we restrict oiirselves to 
exact generahzcd de Sitter solutions in tlie following discussion. We have found four classes of generalized 
de Sitter solutions for (Tp = a-g = and CTp 7^ 0, 0-^ = in Type IIA superstring while one solution for 
(Tp ^ 0, aq = in Type JIB superstring. There is no solution in other cases. 

6.1 ap — aq — 

In this case, we have the same Eqs. (3.1) and (3.2) with 5 = 0. Exchanging 04, 7 and an, 711 and setting q = Q, 
we may solve these equations numerically. In Fig. 6, we give numerical solutions llEi+{au, {i = 1,2) 

with m > with respect to au, and summarize their properties in Table 3. We have two generalized de Sitter 
solutions IIElj- and IIE2± for type IIA superstring, whereas no solution for type IIB superstring, i.e. an = 0. 

6.2 (Tp = 0, Uq ^ (or cTp 7^ 0, cr^ = ) 

It is easy to see that there is no exact solution unless = in which case we have constant Ap = X = 
H^, Aq = aq. Our basic equations reduce to Eqs. (3.15) and (3.16) with (5 = 0. Exchanging 04, 7 and an, 711 
and setting p = 3, q = 6, we solve these equations numerically and find that there is no solutions in the range 
< an < 1/(3 • 2^). Thus we have no solution of this class for both type IIA and type IIB superstring. 

For the case of Up 7^ and ct, = 0, exchanging /i, p and v, q, we find two solutions IIE3± and IIE4-I-. In 
Fig. 7, we give numerical solutions IIEi+(Q!ii, /Xj, fj) {i = 3,4) with fii > with respect to an, and summarize 
their properties in Table 3. 



Table 3: Generalized de Sitter Solutions IIE1+ - IIE4+ with /Xi(Fi) > for < 5s < +00. Five eigenmodes for linear 

pcrturliatious arc also hIiowii. (;r(s,;ni) lueaus tliat tliorc arc in Htalilc modes and /) uustal)le modes. 



Solution 


Property 


Range 


Stability 


3/ii + 7f i 


IIE1+ 
IIE2+ 


1^1 < < fii 
1^2 < < IJ-2 


-0.0047 24 < Q4 < 3~"2-^ (0.5506 < gs < +00) 
-0.0049 49 < 04 < 0.009 447 (0.5751 < Qs < 1-887) 


(ls,2u) 
(0s,3u) 




IIE3+ 
IIE4+ 


At3 = 0, a-p(3) < < z/3 

/LI4 = 0, 5-p(4) < < Z/4 


< 0:4 < 3"'2-'' (0 < ffs < +00) 
< 04 < 3"^2-^ (0 < ffs < +00) 


(2s,2u) 
(3s,lu) 


+ 
+ 



We have two generalized de Sitter solutions IIE3± and IIE4± for type IIA superstring, whereas one solution 
IIE3± for type IIB superstring. 

7 Inflationary Scenario in M-theory 

We are ready for discussion about inflation in M-theory. Since our world is 4-dimensionaL we should also ana- 
lyze our solutions in 4-dimensional Einstein frame, in which the Newtonian gravitational constant is constant. 

7.1 Description of solutions in the Einstein frame 

We have found generalized de Sitter solutions 

a = aoe'^* , b = boe"* , for e = , (7.1) 

and power-law solutions 

a = aoT^ , b = bar" , for e = 1 . (7.2) 

These solutions give the power-law expansion and the exponential expansion in the Einstein frame of our 
4-dimensional spacetime [23] . 
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• Generalized de Sitter solutions (e = 0) 
For u ^0, we have the power-law expansion in the Einstein frame. 



ttE = exp 



Ui 



7U2 



{p - i)m 

qv 



tE = tE^°^ exp 



p-1 



-vt 



(7.3) 



where for > 0, we have ^e^"^ > and t € (—00, 00) is transformed into ^e G (0, 00), whereas for < 
wc have tE*^°'' < and t G (—00, 00) is transformed into t-E, S (—00, 0). For = 0, we have the exponential 
expansion in the Einstein frame 



ttE oc exp[/xfE] , 



tE = t 



(7.4) 



Power-law solutions (e = 1) 

For V ^ — (p — l)/q, we have the power-law expansion in the Einstein frame 



ttE oc tE 



\ = 



{p - i)m + 

(p - 1) + qv 



tE = tE^°^ exp 



p-1 



V t 



(7.5) 



where for v > —(j)—\)/q, we have fE^°^ > and t € (—00, cxd) is transformed into Ie G (0, 00), whereas for 
V < —{p—l)/q, we have tE^^^ < and t G (—00, 00) is transformed into tE G (—00, 0). For v = —{p—l)/q, 
we have the exponential expansion in the Einstein frame 



ue oc exp[(/i - l)tE] 



tE = t , 



(7.6) 



We list the behavior of a scale factor and show the condition for inflation in the Einstein frame in Table 4. 
Note that the values of fi and v in generalized de Sitter solutions (7.1) depend on the choice of the unit. In 
the M-theory, we use the unit of I7I = 1, i.e. mn = 6-i/2(47r)-5/9 _ 0.1818176. If we set TOii = 1, the values 
of n and v in the following tables should be multiphed by the factor 6^/^(47r)^/^ ~ 5.5. On the other hand, 
the power exponent fj, and v in the power-law solutions (7.2) or A in Eqs. (7.3) and (7.5) are dimensionless 
and they do not depend on the choice of the unit. 



Table 4: Behavior of solutions in the Einstein frame; "Condition" means condition for causing inflation in the Einstein 
frame, while the super-inflationary solution behaves as |tE|~''^' for tE — »■ 0-. 











Scale Factor 


Condition 


Range of tE 


Type of inflation 


e = 


iy>0 
j/ = 

u <0 






oe oc explAits] 

flE oc tE"^ 


> 

/X > 
A < 


(0, 00) 
(—00, 00) 
(-00,0) 


power-law 
exponential 
supcrinflation 


e= 1 


v> - 

V = — 

V <- 


{p- 
ip- 
ip- 


1)/? 
1)/? 
1)/? 


Oe oc Ie'^ 

Oe oc exp[(/u — l)tE] 
Oe oc tE^ 


At > 1 

(X> 1 

A < 


(0,00) 
(—00, 00) 
(-00,0) 


power-law 
exponential 
superinflation 



7.2 Conditions for successful inflation and some preferable solutions 

Before we apply our solutions to cosmology, we have to specify what kind of solutions we need. We list 
necessary conditions for successful inflation in our model below. 

(1) iJ.> V and > 0: 

Our four-dimensional universe makes sense only if it is much larger than the internal space, so the external 
space should expand faster than the internal space. Its expansion may not be necessarily inflationary, 
but at least the external space must be expanding in the whole space. Some solutions give an inflation 
in the Einstein frame but the external space shrinks in the original higher dimensions. Such solutions 
are not suitable for a good cosmological model. 
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(2) 60 e-foldings of inflationary expansion; 

We need at least 60 e-foldings of inflationary expansion in the Einstein frame. This may give some 
constraint on the power exponent for a power-law inflation, that is the power exponent should be signif- 
icantly larger than unity. Speciflcally, for generalized de Sitter solution and power-law solution, we find 
that the number of e-foldings A^e in the Einstein frame is given by 



where Uni is the beginning time of inflation in the original frame whereas tan is the ending time. 

(3) (semi)stability against the linear perturbation: 

Solutions which have many unstable modes may not be generic. The most preferable solutions in this 
model are those with only stable modes, which means that the solution is stable against linear perturba- 
tions. However, this kind of solution predicts that inflation never ends because it is stable. For a realistic 
cosmological model, the solution must have small instability for the inflation to end. On the other hand, 
we also want such a solution to be rather generic which requires some sort of stability. This would be 
achieved if the solution contains only one small unstable mode and many other stable modes, and then 
the generic spacetime may first approach this solution and gradually leave it. 

We summarize our solutions in the Einstein frame as well as those in the original frame in Appendix D. 
Using the tables in Appendix D, we shall pick up the preferable solutions for inflation. 

For (5 = 0, although we find exponential expansion of the external space in the original frame, this gives 
non-inflationary power-law expansion in the Einstein frame as described above. There are power-law inflations 
(MPl, MPll, MP12 and MP13) in the past regime. However, the power exponent of them are 1.3 - 2.3, which 
may be too small to solve the flatness and horizon problems, because we do not expect the expansion in these 
solutions continues so long. Thus, these solutions are excluded by the condition (2). 

For S = -0.001, we have six candidates [ME2+, ME3+ (MFl), ME6+ (MF2), MP6 and MP7]. Among 
these, the condition (1) exclude solutions ME3+ (MFl) and MP7 for the internal space expands at the same 
rate as the external space. Solutions ME2+ and MP6 give super-inflation in the Einstein frame. Especially, 
the solution ME2+ has only one unstable mode, and fulfills the condition (3). In this case, we come close to 
the singularity at = 0, but we hope that stringy effects renders the singularity harmless when the curvature 
becomes large. For the remaining solution ME6+, we find an exponential expansion of the external space both 
in the original and the Einstein frames, and the internal space is static, viz. modulus is fixed. From the Table 
2, we also find that this solution fulfills the condition (3) and exists for wide range of S. The radius of the 
external space is arbitrary whereas that of the internal space is bo = 1.893 = 0.3442 tou"^. 

The above solution ME6+ has one unstable mode so this is the most interesting cosmological solutions in 
our criterion. But the instability a^^^ is the same order of magnitude as other eigenvalues of stable modes as 
seen from Eq. (5.16) and appears a little too large to give enough expansion. Specifically, the time scale in 
which this unstable mode becomes important is evaluated as t^s ^ {'^5'^)^^ ^ 1.020. The number of e-foldings 
will be given by A'^e = M^us ~ 0.7905, except for some fine-tuned initial conditions. 

If the eigenvalue of the unstable mode is much smaller than those of other four stable modes, however, a 
preferable solution is naturally obtained for a wide range of initial conditions. We may not need a fine-tuning. 
Although we do not find any value of 6 which gives enough small eigenvalue of an unstable mode, we note 
that our starting Lagrangian has some ambiguity, that is, the forth-order correction term Sw is fixed up to 
the Ricci curvature tensors. We have included additional quartic Ricci scalar term (2.7) in order to take this 
ambiguity into account, but we may still have another kind of correction term including the Ricci curvature 
tensors. To further improve the solutions, we have the possibility of finding more interesting solutions with 
these appropriate extra corrections. 

Another interesting possibility is the following. We have seen that many of our exact solutions have unstable 
modes as well as stable modes. But it is not immediately obvious what happens to any inflationary solutions 
after the solutions decay into other solutions. Our above analysis appears to indicate that we cannot obtain 
big enough e-foldings before the exact solutions decay. However there is a possibility that we may obtain 




(7.7) 
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enough e-foldings after the decay if we can foUow the evolution of our solutions. Indeed we will show in the 
next section that we obtain some numerical solutions which first approach ME64- and give enough e-foldings 
if we fine-tune the initial data. 

7.3 Numerical analysis for generic initial conditions 

To study whether we obtain a sufficient e-foldings in the present model, we have performed the numerical 
calculation around the solution ME6+ for uq = Q, S = —0.1 and cr^ 7^ 0. In this case, Eqs. (2.18) and 
(2.19) give evolution equations for u'l and U2, and Eq. (2.17) gives a constraint equation for seven variables 
w 1, M2, ill, U2, uii U2 and U2, which we used to check the numerical error. Thus, we can give six independent 
initial values and the remainder is given through the constraint equation. Under separate five sets of initial 
conditions, we have performed the numerical calculations for seven dynamical variables. 

In Figs. 8 - 11, we depict five numerical solutions MNl - MN5 whose initial values lie in the vicinity of 
ME6-)-. In Figs. 8 and 9, we show the behavior of solutions in the U1-U2 plane. Generalized de Sitter solutions 
are expressed as a point {ui,U2) = (m? '^) t^iis plane. In these figures, we have the exact solution ME6+ in 
Eq. (4.6) and the future asymptotic solution MFl in Eq. (4.20) (MFl is exact for aq = 0). Every solution 
approaches the ME6+ in the early phase whereas, in the last phase, solutions MNl - MN4 approach the future 
asymptotic solution MFl and the solution MN5 goes to the singularity ui,U2 —00 for finite lengths of time 
as shown in Fig. 9. 




Figure 8: Numerical results around ME6+ in the ui- Figure 9: Numerical results around ME6+ in the U1-U2 
U2 plane for 5 = —0.1. Solutions ME6+ and MFl are plane. The solution MN5 runs away to Ui,U2 —00 in 
indicated by a point in a circle. the last phase. 

We give the initial values as follows: b = 3.7, Ui = 0.28, ua = 0.27, iii = -0.060, U2 = -0.46, u'l = 0.47 for MNl, 

b = 4.4, ill = 0.056, U2 = 0.063, ili = 0.24, U2 = 0.068, ui = -0.15 for MN2, same as MNl except for ui = 0.30 for 

MN3, same as MN2 except for 112 = 0.083 for MN4, b = 5.0, iii = 1.0, U2 = -0.13, ili = -1.6, U2 = 0.20, ui = 4.7 
for MN5. 



In Figs 10 and 11, we show the behavior of the scale factors in the original frame (a, b) and in the Einstein 
frame (oe, 6e), respectively. Here we set qq — bo and t-^it = 0) = 1. In the early phase, every solution gives 
exponential expansion oe c>c e'^-^^^^*'= arising out of ME6+ both in the original frame and Einstein frame, while 
solutions MNl - MN4 give power-law expansion ce oc ^e^'^*^ arising out of MFl in the Einstein frame in the 
late phase. The solution MN4 is not shown in the figures since it does not lead to interesting result. 
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Figure 10: Numerical results in the original frame for 
5 = —0.1. MNi(a) means the plot of In a of the solution 
MNi. We omit the solution MN4. 
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Figure 11: Numerical results in the Einstein frame for 
(5 = —0.1. MNj(a) means the plot of InaE of the solution 
MNi. We omit the solution MN4. 



Despite the analysis in § 7.2, we have enough e-foldings around the solution ME6-|_ with appropriate 
initial values. Especially, numerical solutions MNl and MN2 have interesting features. In the early phase, 
they give the exponential expansion in the external space, and the power-law expansion in the late phase as 
shown in Fig. 11. They do not give enough e-foldings around the solution ME6+ (A^e ~ 30). However the 
solution MFl, which is a stable solution, does also give inflation. Then after leaving the solution ME6+, a 
spacetime approaches to the solution MFl in the late phase, which totally gives 60 e-foldings [N^ = 60) at 
In^E = 35.50 (t = 119.7) for MNl and IniE = 36.62 {t = 119.2) for MN2. The scale of internal space becomes 
Ra = b = 2.182 x 10* = 3967 tou-i for MNl and Rq = 2.995 x 10"* = 5446 jtiii^i for MN2 both in the original 
and Einstein frames. After inflation, if the internal space settles down to static one, the present radius of extra 
dimensions is given by Rq. Since this is slightly larger than the fundamental scale length, we may adopt the 
model of large extra dimensions, which was first proposed as a brane world by Arkani-Hamed et al. [25]. In 
this model, the 4D Planck mass is given by TO4^ = i?o^mii^. We then flnd 

TOii = 2.543 X 10"^^ m4 = 619.5 TeV , for MNl , (7.8) 

mil = 8.392 X 10"^"* m4 204.4 TeV , for MN2 . (7.9) 

This is our fundamental energy scale. The scale of extra dimensions is Rq = 3967 mn^^ — 6.403 TeV"^ for 
MNl and Rq = 5446 mn-^ = 26.64 TeV"^ for MN2. 

Although numerical solutions MNl and MN2 give interesting scenario, their behaviors highly depend 
on choice of initial values. For example, numerical solutions MN3 and MN4 have almost the same initial 
conditions as those of MNl and MN2 (see Fig. 8), but their final behaviors are very different (Fig. 11). 
We can find 60 e-foldings at IntE = 49.97 {t = 83.47) for MN3, and the scale of internal space becomes 
Rq — 1.349 X 10^ = 2.463 x 10'"' nin^^. Thus, if the extra dimension is so large, the IID Planck mass turns 
out to be TOii = 0.3285 GeV, which is excluded by high-energy experiments. 
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8 Concluding remarks 



In this paper, we have examined cosmological sohitions of the effective theories of M-theory and superstrings 
with special attention to the de Sitter-like and power-law expansions. In order to evade the no-go theorem in 
this setting, wo have included higher order quantum corrections. Wc have also taken the ambiguity arising 
from the field redefinition into account. We have found that there is an interesting solution ME6+ in which 
the external space is expanding whereas the size of the internal space is static. This is interesting in another 
respect that this may be regarded as "moduli stabilization" by the higher order corrections. It is true that 
this solution does not realize the usual "moduli stabilization" where all moduli are fixed and stable because we 
have one unstable mode in the fluctuations around this solution. But if the solution were completely stable, 
the inflation would not end. On the other hand, it is necessary that the size of the internal space does not 
grow much during inflation. So the feature of this solution is desirable for giving successful inflation. Naively 
it may be expected that there are many inflationary solutions. It is somewhat surprising to find that we have 
found few solutions with desirable features. 

Quite interestingly, we find some numerical solutions which first approach an exact infiationary solution 
and then give enough e-foldings if we fine-tune the initial data. The key to understand this result is that 
we cannot see how the solution evolves after the exact solution decays if we are looking only at the analytic 
solutions. We were able to follow the evolution by the numerical analysis. Even though it appears that this 
may need some special initial conditions, we find that this opens another possibility of achieving successful 
inflation and is a very interesting result. 

There is also still ambiguity in other terms involving Ricci tensors and scalar curvatures in the effective 
theory, which we have not analyzed. So rather than taking our results literally, we should understand these 
results as an indication that such inflationary solutions are possible in this direction. It is certainly an inter- 
esting problem to see if our results may be further improved by including other possible quantum corrections. 
We also have to look at the effects of fluxes, which play important roles in inflationary models with an effective 
potential [26] . These may lead to successful inflation insensitive to initial conditions. 
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A Explicit form of actions 

We use the following notation throughout this paper. 

{p — m)n ^ {p — m){p — m — l){p — m — 2) ■ ■ ■ (p — n), 
{q-m)n = {q-m){q-m-l){q-m-2)---{q-n), 

X = ill — uoiii + ill, Y = U2 — U0U2 + 

Inputting our ansatz (2.15) into actions (2.2)-(2.9), we can write down the Lagrangians in terms of uq, ui 
and W2- The explicit forms of actions are listed here, where actions (A. 2) - (A. 5) are integrated by parts. 

(1) Einstein-Hilbert action (n = 1) 

E2 = e-^'"' [piAp + qiAg - 2{piUi^ + pquiU2 + 51^2^)] • (A.2) 

(2) Gauss-Bonnet action (n = 2) 

£14 = e-^"° [p3^p^ + 2piqiAgAp + q^Ag^ 

- 4Ap(p3iii^ +p2quiU2 +PiqiU2^) - AAgipiqiiii^ +pq2UiU2 + 93^2^) 

+ ^(2p3Mi* + 2p2qui^U2 + 3pigiWi^U2^ + 2pq2U\U2^ + 2g'3«2*)] . (A.3) 

(3) Lovelock action (n = 3, 4) 

= e^^"" [p5Ag^ + SpsqiAp'^Ag + 3^193^^71,2 _^ ^^^^3 _ QAp'^ip^ui'^ +P4gwiW2 +P3?iM2^) 

- 6Aq2(pi53Mi^ + + g5W2^) - 12ApAg{p3qiUi'^ + p2q2UiU2 + PiqsU2^) 
+ AAp{2p^Ui^ + 2p4qui^U2 + 3p3qiUi^U2^ + 2p2'72MiW2^ + 2^193112'') 

+ 4Aq(2p3(?i?j,i"' + 2p2'?2Ui^U2 + 3piq3Ui'^U2^ + 2pq4UiU2^ + 2q:^U2^) 

+ g(2p5Wl*^ + 2piqUi^U2 + 5p3q'iWi^U2^ + ^P2q2U\'U2^ + 5piQ3'«l^W2'^ + 2pq^UiU2^ + 2g5U2^)] , (A.4) 

£■8 = e"^"" [pt^p" + ^pmiAp^Ag + Qp^q^Ap^Ag^ + 4^195^^^,^ + q-jAg^ 

- SAp^{p-jUi' +pQquiU2 +P59iW2^) - SAg^ {piqr,ui' +pqeuiU2 + 97^2^) 

- 24:Ap'^Ag{p5qiUi'^ + p4q2UiU2 + P3q3U2'^) - 24:ApAg'^ {psq^iii'^ + p2q4UiU2 + Piq5U2'^) 
+ 8Ap'^{2p^Ui^ + 2pequi^U2 + 'ip^qiUi'U2^ + 2p4q2UiU2^ + 2p3g3U2^) 

+ 8Aq2 (2^393^1'' + 2p2qAUi^U2 + ipiqbUi^U2^ + 2pqQUiU2^ + 2q-jU2'^) 

+ l&ApAg{2p^qiUi' + 2p4q2Ui^U2 + 'ip3q3Ui'U2^ + 2p294UiU2^ + '^PiqbU2^) 
32 

- —Ap{2p'jUi^ + 2pQqUi^U2 + ^PhQlUl^U2^ + 5p4g2Wl^W2^ + ^P3q3Ul^U2'^ + 2p254WlW2^ + 2piq^U2^) 


32 

- -^Ag{2p^qiui^ + 2p4q2Ui^U2 + 5p3g3Ui*U2^ + 5p294Wi^W2^ + 5p5q3Ui'^U2'^ + 2pqQUiU2^ + 2qrU2^) 
16 

+ -^{S>P7Ui^ + Spequi"^ U2 + 28p5(7iUi^U2^ + 28p4q2Ui'^ U2^ + 35p3q3Ui'^U2'^ 
35 

+ 28p2g4Mi^U2^ + 28piq5Ui^U2^ + 8pqeUiU2'^ + 8g7M2*)] • (A.5) 
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(4) Sw action 

Components of the Weyl tensor defined in (2.10) are given by 

C*at5 = -77^ -TTT^ —pBlQab, 



{D- 


■1){D- 


2)^ 








{D 


-1){D 


- 2 




g-2«o 




{D- 


■1){D- 


2)^ 




g-2«o 




{D- 


-1){D- 


2f 




g-2uo 


1 



bed 

where i, . . . run over external space, and a,b,... over internal space, respectively, and we have defined 

Bi = {D- 3)(X -Y)-{;p- l)Ap + {q - + [p - q)uiU2 , 
B2 = -2{X -Y) + {q + l)Ap + {q- l)Ag - 2quiU2 , 
S3 = 2{X -Y) + {p- l)Ap + {p+ l)Aq - 2puiU2 , 

B4 = ip- q){X -Y)-{p- l)qAp - p{q - l)Ag + {2pq -p- q)uiU2 . (A.7) 

Although wc have four quantities Bi, B2, B^ and i?4, only two of them are independent. Actually, the 
tracelessness of the Wcyl tensor gives the following relations. 

Bi + ip- l)B2 + B4 = , -Bi + {q- l)B3 + B4 = . (A.8) 

Taking independent variables as B2 and in order to consider the case p ~ 1 or q = I at the same time, we 
find that the other variables Bi and B4 are given as the following with respect to B2 and B^. 

B^=-^[{p-l)B2-{q-l)Bs] , B4 = -^[{p-l)B2 + {q-l)Bs] (A.9) 

Substituting Eqs. (A.9) and (A. 6) into the action Eq. (2.6), we have the Lagrangian with respect to B2 and 
B3 as follows. 

" 16(jj-i)4(j^-2)4 - An2pgB2^Bs + 2n3pgB2^-B3^ - 4n2gpB2Bs^ + nigpSg^] (A.IO) 

where we define riipq, n2pq and n^pq as 

nipq = (j? - 1) {-{p - \)\p - 2) + (p - l)3(p2 + 3p _ 7)g 

+ 2{p - l)2(p2 _ 4p + 7) + (p3 _ ;^3p2 ^ _ ^^^^Sj ^ 

n2pq = (P - l)'(g - 1) [-3p^ + 3/ + p2 + p - 2 + (p - l)(p3 + 2p2 + 5)g 

+ 2(p(p - 1)2 - 4)92 + (p _ 3)(p - 5)g3] , (A.12) 

nzpq = {p- l){q - 1) [p{p^(3g2 - 18(z + 23) - ^^(Qg^ - 16g + 13) - p(llp + 1) - 3} 
+ q{qH^P'^ - 18p + 23) - q'^{9p^ - 16p+ 13) - g(llg+ 1) - 3} 

+ 2{3(pqf + 5ipqf + 8pq - 3}] , (A.13) 
while niqp and n2qp are given by changing p for q in Eqs. (A. 11) and (A.12), respectively. 

(5) i?* action 

Lr^ = e-7"o+P«i+9«2 ^2pX + 2qY + piAp + q^Aq + 2pquiU2] ^ . (A.14) 
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B Field equations 



Taking variation of the actions (A.2)-(A.14), we find the basic equations (2.17) - (2.19), where each term is 
summarized here according to which action it originates from. The expUcit forms of each term in the field 
equations are listed hero: 

(1) EH action (n = 1) term 

Fi = aie""° [piAp + qiAg + 2pquiU2] , (B.l) 

/(P) = aie-"« [{p - l)2Ap + giA, + 2{p - l)quiU2] , (B.2) 

= aie-"° [piAp + {q- 1)2^, + 2p{q - 1)^1^2] , (B.3) 

g'f^ = 2{p - l)aie-"° , g[''^ = 2{q - 1)016""° , (B.4) 

h^f'^ = 29016-"" , h^^^ = 2^016-"" . (B.5) 

(2) Lovelock action (n = 4) term 

F4 = 046-^"" [p7^p + 4^591^3^, + epsqsAlAl + 4piq5ApA^g + qjA\ 

+SuiU2{p&qAl + ipiq2ApAq + ip2qiApA^q + pq^A^^) 

+24ulul{p5qiAl + 2psq3ApAq + piq^A^) + 32wfw^(p4g2^p +P2«4^g) + IQpsqsutut] , 



/i^) = a4e-'"« [{p - 1)8^^ + 4(p - l)eqiAlAq + 6{p - lUqaA^A^ + 4(p - l)2q5ApAl + q^A^ 

2 



(B.6) 



+8uiii2 {{p - l)7qAl + 3(p - l)5q2AlAq + 3(p - l)3,qiApAl + {p- l)g6>i^) 
+24^2^2 ((p _ 1)6^1^2 + 2{p - l)4q3ApAq + {p- 1)295^9) 

+32ulul {{p - l)592Ap + {p- 1)394^9) + 16(p - l)4q3utut] , (B.7) 
/i'^ = 046-^"° [{q - 1)8^^ + 4(9 - l)ePiA^qAp + 6{q - l)4P3^p^^ + 4(9 - 1)2^5 + Pr^t 
+8uiU2 {{q - IhpAl + 3{q - 1)5^2^^^^ + 3{q - 1)3^4^,^^ + {q - l)P6^p) 

+24^2^2 ((g _ l)epi^2 ^ 2{q - l)4P3^p^g + (? - l)2P5^p) 

+32u?ui ((g - 1)5P2^9 + {q- 1)3^4^^) + 16(9 - l)4P3wfw|] , (B.8) 
ffi^^ = 8(p - 1)046-^"° [(p- 2)7^3 +3(^-2)591^2^, + 3(^-2)393^^^2 + 95^3 

+6uiU2 ((p - 2)69^2 + 2{p - 2)4q2ApAq + {p- 2)94^^) 

+ 12«2u2 {{p - 2)5<ZlAp + {p- 2)3(?3^g) + 8(p - 2)4q2ulul] , (B.9) 

ffi'^ = 8iq - l)a4e-'"« [(g - 2)7^^ + 3(g - 2)5^1 A^A^ + 3(9 - 2)3^3^^ + P^^l 
+6wiU2 ((g - 2)6P^2 _^ 2(9 - 2)4P2ApAg + (q- 2)p4Al) 

+I2ulvl {{q - 2)5PiAq + {q- 2)apaAp) + 8{q - 2)4P2M-i?i2] , (B-IO) 
h'f^ = 8<za4e-^"" [{p - l)eAl + 3(p - 1)4(9 - 1)2^^^, + 3(p - 1)2(9 - l)iApAl + (9 - 1)6^^ 
+6uiU2 {{P - 1)5(9 - l)Al + 2{p - 1)3(9 - l)3ApAq + (p - 1)(9 - 1)5^^) 
+12ulul {{p - 1)4(9 - l)2Ap + ip~ 1)2(9 - 1)4^,) + 8(p - 1)3(9 - Ihulul] , (B.ll) 
4'^ = 8pa4e-^"° [(9 - 1)6^9 + 3(9 - 1)4(P - 1)2^2^1^ + 3(9 - 1)2(P - 1)4^,4 + {p- l)eAl 
+6«iU2 ((9 - 1)5(P - l)Al + 2(p - 1)3(9 - l)3ApAq + (9 - l)(p - 1)5^) 
+12ulul ((9 - l)4(p - 1)2 A + (9 - 1)2(P - l)4Ap) + 8(p - 1)3(9 - l)3ulul] . (B.12) 
(3) Sw action term 



V = ^e-f"i-«"^ 



dB2 ' ' 9^3 
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dLw dL 



■w 



.W-2,,e'<..-«..{(, + l,^ + (p-l,^} 

\ dLw / . \ dL 



)}] , (B.13) 



„d ( / , , \ dLw ( . , , . . \ dLw \ 

_ rf^ f dLw _ dLw \ 
dt'^ I dB2 dBsii ' 



(B.14) 



[qL 



w — 2aqe 



2(«o-'i 



„ d f / . . , , X . \ dLw ( . . , 1 N . \ dLw ~\ 



^d^(dL 



w 



dL 



w 



where 



dLw pq e 



-7ua+pui+qu2 



dB2 

dLw pq e 



A{D ~lf{D -2Y 

7«0+PMl+g«2 



(4) i?'' action term 



pipg-B2 — 3n2pq-B2 i?3 + n3pgi?2-B3 ~ n2qpBz \ , 
[—'n2pqB2^ + n3pqB2^ B3 — 3n2gp-B2-B3^ + niqp-B3^] . 



(B.15) 

(B.16) 
(B.17) 



■ iJ4 



where 



pLn^ - 2ape2("o-»i) 



91/ /f 4 d 
dAp 



'9 



(B.18) 



dLjii . 9-L/J4 9ii?4 \ f dL 



dX 



- 2ui 



dAr^ 



dill 



dt^ 



, . „ . , dLni „ . dLni dLjiA. \ d 

[Uq - 2U2)—^ ^U2-^, — I + — 



dY 



dA„ 



dU2 



a 
d^ 



dX 
(B.19) 

dLRi 



dY 
(B.20) 



dX 
dLijA 

dAp 

dLni 

dill 



-7«0+P«l+9M2 



Spw 



dL 



dY 
dA„ 



-7uQ+pui+qU2 



8qR^ 



-7«0+P«l+9M2 



8pqu2R 



dU2 



R = 2pX + 2qY + piAp + qiAq + 2pquiU2 



(B.21) 



(B.22) 



C Inputting our ansatz into solutions 

In order to find solutions, we assume 

Uq = et , ui = fj,t + In oq , U2 = i^t + In 60 



(C.l) 
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Inserting this form into the above equations (Eqs. (B.l) - (B.15)) and setting 



A* +ape 



2(e-/i)t 



,2 ' 



?)2 ' 
"0 



X = fj,{n-e) , Y = v{v-e) , 
we obtain the following explicit equations: 
(1) EH action (n — 1) term 



(p) 



aie [piAp + qiAq + 2pq(iv\ 
aie~'* [{p - l)2Ap + q^Aq 



n 



(9) 



9i 



(p) 



Ap) 



(2) Lovelock action (n = 4) term 



aie-"* [piAp + 



2(p - 1)?^^ , 
{q-l)2Ag + 2piq-l)fiiy], 
2{q-l)aie~'\ 



h[^'^ = 2paie-'\ 



Fi = a^e-^'* [prA^p + ipsQiAlAg + 6p3«3^p^^ + "^PiqsApA^ + q^A] 



+8fii^{peqAl + 3p4q2AlAq + 3^294^^^, +OT6^q) + 24:fi^i^^{p5qiAl 
+2p3q3ApAq + piq^Al) + 32^i^iy^{p4q2Ap + p2q4Ag) + 16^393^^'^''] , 



ip) 



-7et 



{p - 1)8^ + 4(p - 1)691 A^A, + 6(p - 1)4934^2 + 4(p - l)2q5ApA: 
+q'jA\ + i^iv{{p - l)7qAl + 3(p - 1)592^2 A, + 3(p - l)3qiApAl 

+ {p - 1)^6^9} + 24m'z^' {(P - l)6ai4 + 2b - l)4«3ApA9 + (P - 1)2«5A^} 

+32m^i/3 {(p - 1)5^2^? + {p- 1)394AJ + 16(p - 1)493^^^^*] , 

8(p - l)a4e-^^* [{p - 2)jAl + 3(p - 2)5^1 A^A, + 3(p - 2)393^^^^ + q^A^ 

+&liy{{p - 2)QqAl + 2{p - 2)iq2ApAg + {p - 2)94^9} 

+12^V2 {{p - 2)5^1 Ap + {p- 2)393^9} + 8(p - 2)492^^2^^] , 

Sqa^e-^'* [{p - l)eAl + 3{p - 1)4(9 - 1)2^^^, + 3(p - 1)2(9 - ^UApA^ 

+ iq " 1)6^^ + 6^2^ {(p - 1)5(9 - l)Ap + 2(p - 1)3(9 - l)3ApAg 

+{P - 1)(9 - 1)5^2} + 12fi^u^ {(p - 1)4(9 - l)2Ap + {p- 1)2(9 - 1)4^} 

3, ,3 



+8(p- 1)3(9-1)3/^^1'' 



Aq) 



= 0:46 



-7et 



(9 - 1)8^9 + 4(9 - l)&PiAlAp + 6(9 - l)iPiAlAl + 4(9 - 1)2^5^,4 
+pjAl + {(9 - l)7Mg + 3(9 - l^PiA^Ap + 3(9 - 1)3^4 A A^ 

+ (9 - 1)P64} + 2^2^^ {(9 - l)6PlA' + 2(9 - 1)4^3^^^, + (9 - l)2P5Ap} 
+32/zV3 {(g - 1)5^2^9 + (9 - l)3P4Ap} + 16(9 - 1)4^3^*1^1 , 

8(9 - 1)^46-'^* [(9 - 2)7A3 + 3(9 - 2)5Pi A^Ap + 3(9 - 2)3P3AgAl+p5Al 
+6^lly{{q - 2)epAl + 2{q - 2)4^2 A, + (q - 2)^4^^} 

+ 12/^2^2 _ 2)5PlAg + {q- 2)3P3Ap} + 8(9 - 2)iP2H^iy'^] , 



-7€t 



(9 - 1)6^^ + 3(9 - 1)4(P - l)2A2Ap + 3(9 - 1)2(P - l)4AgAl 



+{p - l)eAl + 6fip{{q - lUp - 1)A2 + 2{p - 1)3(9 - l)3ApAg 
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+ iq - l){p - 1)5^2} + 12Ai2z.2 {(q _ _ i)^Ag + (g - IUp - l)4Ap} 
+8(p-l)3((Z-l)3/x'i/3], (C.ll) 

(3) Sw action term 

= 16(D^- 1)4 (J - 2)4 [-^^^ - - ^)(-^^ + '?'^)^- 
+ 2ape2(^-^)*{iV+ + N_ + {e- ^i){{q + 1)7V_,2 + (p - l)iV-,3) } 

+ 2age2(^-'')*{iV+-iV_ + (e-i.)((g-l)iV_,2 + (p+l)iV_,3)}] • (C.12) 
^^^"^ = 16(D-\)4(Z)-2)4 [^^^ - + + 9i^)((-6e + (p - + qu)N_ + - i/)7V+) 

- 2CTpe2(^-'')*{Af+ + iV_ + 2(e - + (2p - 3)m + 2qu){{q + 1)7V_,2 + (p - l)A^-,3) 
+ (M - + l)iV+,2 + (P - l)A^+,3))} 

- 4(e - !/)a,e2(^-'')*{(-lle + (2p - 1)^ + 2{q - l)u){{q - 1)N_^2 + (p+ l)A^-,3) 
+ ,y){{q ~ l)7V+,2 + {p + l)iV+,3)} 

- 8(e - ^ira/e^^^-^^^'{{q + lfN_^22 + 2{p - l){q + 1)7V_,23 + {p - l)'iV_,33} 

- 16(e - M)(e - z/)apage2(2^-^-'^)*{(g2 _ i)iv_,22 + 2(p9 + 1)7V_,23 + (p + l)'iV_,33} 

- 8(e - i.)2a,2e4(--)*{(g - 1)27V_,22 + 2{p + l){q - l)iV_,23 + (p - l)'iV_,33}] , (C.13) 

Fw^"^ = i6(dI^i)4(^*_2)4 [^^^ + 2(-7e + p/i + g^.) ((-6e + + (? - l)u)N. + (/x - i/)7V+) 

+ 4(e - ;u)ape2(^-'')*{(-lle + 2(p - + {2q - l)u){{q + 1)7V_,2 + (p - l)7V-,3) 
+ (M - u){{q + l)7V+,2 + (p - l)Ar+,3)} 

- 2age2(^-")*{Ar+ - Af_ - 2(e - !/)((-lle + 2p/i + (2g - 3)u){{q - 1)7V_,2 + (p + l)A^-,3) 
+ iy){{q - l)7V+,2 + (p + l)iV+,3)) } 

+ 8(e - M)'a/e4(-^)*{(<z + 1)27V_,22 + 2(p - l)(g + l)iV_,23 + (p - l)'A^-,33} 

+ 16(e - Ai)(e - i^)apage2(2e-M--)t|(^2 _ ^ 2(p(7 + l)Af_,23 + (p + l)'iV_,33} 

+ 8(e - !/)2a,2e4(^--)*{(g - 1)27V_,22 + 2(p + l){q - 1)7V_,23 + (p - l)'^-,33}] • (C.14) 

where 

Lw = nipgB2^ - An2pqB2^B3 + 2n3p,B2^B3^ - AniqpBiBs^ + niqpB3^ , (C.15) 

Ma = 4 [nipqSa^ - 3n2pgB2^S3 + nspgB2Bs'^ - naqpBs^] , (C.16) 

M3 = 4 [-n2pqB2^ + nspqB2^Bs - 3n2qpB2Bs^ + m^pBs^] , (C.17) 

7V+ = qM2 + pMa 

= A[{qnipq -pn2pq)B2^ - {3qn2pq - pn3pq)B2'^ B3 

+ {qn^pq - 3pn2qp) -82-83^ - (gn2gp - pnigp)S3^] , (C.18) 

Ar_ = M2 - M3 

= 4[(nipq + n2pq)B2^ - (3n2pg + n3pq)-B2^-B3 + (nspq + 3n2gp)-B2-B3^ - (n2gp + niqp)B3^] , (C.19) 

7V+,2 = 4[3(gnipq - pn2pq)B2'^ - 2{3qn2pq - pn3pg)-B2-B3 + {qn3pq - 3pn2qp)B3^] , (C.20) 

= 4[-(3gn2pg - pn3pq)B2^ + 2{qn3pq - 3pn2qp)B2B3 - 3{qn2qp - pnipq)B3^] , (C.21) 
7V_,2 = 4[3(nipg + n2pq)B2^ - 2{3n2pq + n3pq)B2B3 + (nsp, + 3n2gp)B3^] , (C.22) 
7V_,3 = 4[-(3n2pg + n3pq)B2^ + 2(n3pg + 3n2qp)B2B3 - 3{n2qp + nigp)B3^] , (C.23) 
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= 8 [3(nip, + n2pq)B2 - (Snzp, + nsp,)^,^] , (C.24) 
^-,23 = 8[-(3n2p9 + nzpq)B2 + (nspq + 'in2qp)Bz] , (C.25) 
7V_,33 = 8 [{nspg + 3n2gp)B2 - 3{n2gp + nigp)B3^] • (C.26) 

(4) action term 

Fn^ = Se-'^'^R^ [-IR^ + 8(-7e + p/i + qv){p^i + qv)R + IQdpC^^''-^'^* [pR + 6pi(e - At)(PAt + 

+ 16a,e2(^-'^)*{9i? + 6(?i(e - z/)(pM + qv)]] , (C.27) 

Ffl4(P) = Se-^"R [R^ - 8(-7e + pfi + qi^){6e + fi)R^ 

+ 8{p - l)ape2('-^)*{-i?2 + 6p(e - fi){-13e + (p - 1)^ + qi^)R + 24p{e - 

+ 46qi{e - i/)5-ge2(^-'')*{(-13e + {p - 1)/Lt + qi^)R + 4(e - ly)} + I92pi^{e - nfa/e^''^-''^^R 

+ 384pigi(e - ii){e - t^)ap^qe2(2e-M-^)t^ + 192gi2(e - iyfdg'^e^^'-''^*R\ , (C.28) 

Ffl4(«) = 66-"^'* R [R^ - 8{-7e + pn + gi/)(6e + i^)R^ 

+ 46pi(e - i^)ape^^'-''^*{{-13e + pii + {q - l)iy)R + 4(e - n)} 

+ 8{q - l)CTge2(^-")*{-^2 + 6q{e - i/)(-13e + + {q - l)iy)R + 24g(e - lyf} 

+ 192pi2(e - /x)2aj,2g4(e-^)t^ ^ 384pigi(e - /x)(e - i/)CTp^ge2(2^-''-'''*^ 

+ 192gi2(e-z/)2CTg2e*(^-")*^] , (C.29) 

where 

^ = 2pX + 2qY + piAp + q^Ag + 2pqiiu . (C.30) 



D Summary Tables of solutions in M-theory 

Here we summarize our solutions for 6 — 0, —0.001 and —0.1 in tables. In the last columns of the tables, we 
include the type of two spaces {dSp,dSg). K means the kinetic dominant space, in which the curvature term 
(cTp, or (Tg) is cither zero or can be asymptotically ignored. M denotes the Milne-type space, which is described 
by ds"^ = —dt^ + t^dSp + • • • with cTp = — 1, or ds'^ = —dt^ + • • • + i^ds^ with ag = —1. Similarly, we define 
a constant curvature space C by cTp = 1 or cTg = 1, and So and S± are static spaces with zero curvature and 
positive (or negative) curvature, respectively. 
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Tabic 5: Exact solutions for 5 = 0. e = and 1 correspond to generalized de Sitter solutions (a ^ e^*, 6 e"') and 
power law ones (a ~ t^, b ~ t"), respectively. A is a power exponent of power law solutions in the Einstein frame 
(ob ^ t'^). K, S±, So, and M mean a kinetic dominance, a static space with positive (or negative) curvature, a flat 
static space, and a Milne-type space, respectively. There is no infiationaxy solution in the Einstein frame. 



Solution 


e 


ffp 


0-q 




V Oo 


bo 


A 




Type 


ME1± 











±0.1047 


TO.9367 • • • 




0.9681 


TO.2676 


K K 


ME12 


1 





-1 





1 


1 


0.7778 


0.7778 


So M 


MEi:i 


1 


-i 





i 


i 




i 





M So 



Table 6: Future asymptotic solutions {t — » oo) for 5 = 0. There is no inflationary solution in the Einstein frame. 



Solution 


e 








u 


Oo 


6o 


A 


<^i 


tE 


Typo 


MF6 


1 








0.5583 


-0.0964 






0.3333 


-0.1455 


— > OO 


Kasner 


MF7 


1 








-0.3583 


0.2964 






0.3333 


0.1455 


— > oo 


Kasner 


MF8 


1 


-1 


-1 


1 


1 


0.4714 


0.8165 


1 


0.2222 


— > oo 


M M 



Table 7: Past asymptotic solutions {t -oo) for .5 = 0. MPl, MP4, MP5, MPll, MP12 and MP13 are inflationary 
solutions in the Einstein frame. 



Solution 


t 


Cp 




M 




Oo 


&o 


A 


01 




Type 


MPl 


1 





0, ±1 


1.588 


0.3193 






1.278 


0.1508 


~ 


K K 


MP2 


1 


0, ±1 


0, ±1 


0.7336 


0.082 88 






0.7935 


0.064 25 


~ 


K K 


MP3 


1 


0, ±1 


0, ±1 


0.7225 


-0.1667 






0.3335 


-0.4004 


~ 


K K 


MP4 


1 


0, ±1 


0, ±1 


0.6221 


-0.4004 






1.942 


0.9975 


~ 


K K 


MPS 


1 


0, ±1 


0, ±1 


0.1003 


-1.701 






1.182 


0.3434 


~ 


K K 


MP6 


1 


0, ±1 


0, ±1 


0.022 04 


0.9906 






0.7811 


0.2218 


~ 


K K 


MP7 


1 


0, ±1 


0, ±1 


-0.030 14 


0.6209 






0.6754 


0.1957 


~ 


K K 


MPS 


1 


0, ±1 


0, ±1 


-0.3353 


0.8502 






0.6641 


0.2139 


~ 


K K 


MP9 


1 


0, ±1 


0, ±1 


-0.6685 


0.6343 






0.4818 


0.1970 


~ 


K K 


MPIO 


1 


0, ±1 





-0.9380 


2.573 






0.8063 


0.2571 


~ 


K K 


MPll 


1 





-1 


6.680 


1 






2.262 


0.2222 


~ 


K M 


MP12 


1 





-1 


6.086 


1 






2.1302 


0.2222 


~ 


K M 


MP13 


1 





-1 


2 


1 






1.222 


0.2222 


~ 


K M 


MP14 


1 


0, ±1 


-1 


0.4818 


1 






0.8848 


0.2222 


~ 


K M 


MP15 


1 


0, ±1 


-1 


0.072 89 


1 






0.7940 


0.2222 


~ 


K M 


MP16 


1 


1 


0, ±1 


1 


0.2682 






1 


0.1383 


~ 


M K 


MP17 


1 


1 


0, ±1 


1 


-9.178 






1 


0.2949 


~ 


M K 


MP18 


1 


-1 


0, ±1 


1 


0.9318 






1 


0.2187 


~ 


M K 


MP19 


1 


-1 


0, ±1 


1 


-0.1225 






1 


-0.2144 


~ 


M K 


MP23 


1 


±1 


-1 





1 


1 




0.7778 


0.2222 


~ 


S± M 


MP24 


1 


-1 


±1 


1 







1 


1 





~ 


M S± 


MP25 


1 


1 


-1 


1 


1 


0.6495 


0.8898 


1 


0.2222 


~ 


C M 


MP26 


1 


-1 


-1 


1 


1 


0.8413 


0.091 06 


1 


0.2222 


~ 


M M 


MP27 


1 


-1 


-1 


1 


1 


0.4780 


2.284 


1 


0.2222 


~ 


M M 


MP28 


1 


-1 


-1 


1 


1 


0.4166 


0.6216 


1 


0.2222 


~ 


M M 
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Table 8: Exact solutions for 5 = —0.001. e = and 1 correspond to generalized de Sitter solutions (a ~ e***, b ~ e"*) 
and power law ones {a ^ t*^ , 6 ~ r"), respectively. A is a power exponent of power law solutions in the Einstein frame 
(as ^ tB).K, S±, So and M mean a kinetic dominance, a static space with positive (or negative) curvature, a flat static 
space, and a Milne-type space, respectively. ME3± and ME6+ are inflationary solution in the Einstein frame, but our 
3-space is shrinking in ME3_. 



Solution 


e 


<7p 






V 


ao 


6o 


A 


01 


Type 


ME2± 











±1.438 


TO.06707 






-5.126 


±0.01916 


K K 


ME3± 











±0.4029 


±0.4029 






1.286 


±0.1151 


K K 


ME6± 








1 


±0.7751 







1.893 







K S+ 


ME8± 





1 








±0.4918 


1.0472 




1 


±14.05 


S+ K 


ME9± 





1 








±0.4016 


1.220 




1 


±0.1167 


S+ K 


ME10± 





-1 








±0.4084 


0.7124 




1 


±0.1147 


S- K 


ME12 


1 





-1 





1 




1 


0.7778 


0.2222 


So M 


ME13 


1 


-1 





1 





1 




1 





M So 



Table 9: Future asymptotic solutions (t oo) for 5 = —0.001. MFl and MF2 are inflationary solution in the Einstein 
frame. 



Solution 


e 


CTp 






V 


Oo 


6o 


A 


01 


ts 


Type 


MFl 





±1 


±1 


0.4029 


0.4029 






1.286 


0.1151 


— > 00 


ME3+ 


MF2 





±1 


1 


0.7751 







1.893 







— > oo 


ME6+ 


MF3 





1 


±1 





0.4918 


1.0472 




1 


14.05 


—^ OO 


ME8+ 


MF4 





1 


±1 





0.4016 


1.220 




1 


0.1167 


— »■ oo 


ME9+ 


MF5 





-1 


±1 





0.4084 


0.7124 




1 


0.1147 


— > OO 


ME10+ 


MF6 


1 








0.5583 


-0.0964 






0.3333 


-0.1455 


—^ 00 


Kasner 


MF7 


1 








-0.3583 


0.2964 






0.3333 


0.1455 


— > oo 


Kasner 


MF8 


1 


-1 


-1 


1 


1 


0.4714 


0.8165 


1 


0.2222 


— »■ oo 


M M 
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Table 10: Past asymptotic solutions (t -oo) for 5 = -0.001. MPl, MP2, MP7, MP9, and MP18 are inflationary 
solution in the Einstein frame, although our 3-space is shrinking in MPl and MP2. 



Solution 


e 


(Tp 




A* 


V 


ao 


&o 


A 




Ie 


Type 


MPl 





±1 


±1 


-0.4029 


-0.4029 






1.286 


-0.1151 


— > — CX3 


ME3+ 


MP2 





±1 


1 


-0.7751 







1.893 







— oo 


ME6+ 


MP3 





1 


±1 





-0.4918 


1.0472 




1 


-14.05 


—oo 


ME8+ 


MP4 





1 


±1 





-0.4016 


1.220 




1 


-0.1167 


— » —00 


ME9+ 


MP5 





-1 


±1 





-0.4084 


0.7124 




1 


-0.1147 


— oo 


ME10+ 


MP6 


1 





0,±1 


121.2 


-5.488 






-5.603 


0.3014 


~ 


K K 


MP7 


1 








27.08 


27.08 






1.272 


0.2827 


~0 


K K 


MPS 


1 





0,±1 


26.66 


-37.15 






0.8011 


0.2879 


~ 


K K 


MP9 


1 





0,±1 


2.610 


-0.1187 






3.756 


-0.2032 


~ 


K K 


MPIO 


1 


0,±1 


0,±1 


0.7376 


-0.086 31 






0.6240 


-0.1237 


~ 


K K 


MPll 


1 


0,±1 


0,±1 


0.7268 


-0.1514 






0.4187 


-0.3221 


- 


K K 


MP12 


1 


0,±1 


0,±1 


0.1909 


0.1395 






0.4564 


0.093 75 


~ 


K K 


MP13 


1 


0,±1 


0,±1 


0.1548 


0.1548 






0.4519 


0.1004 


~ 


K K 


MP14 


1 


0,±1 


0,±1 


0.1201 


0.1699 






0.4483 


0.1066 


~ 


K K 


MP 15 


1 


0,±1 


0,±1 


-0.7576 


0.6250 






0.4486 


0.1961 


~ 


K K 


MP16 


1 


0,±1 





-1.162 


1.498 






0.6537 


0.2399 


~ 


K K 


MP17 


1 


0,±1 


0,±1 


-2.408 9 


0.5986 






-0.1009 


0.1934 


~ 


K K 


MP18 


1 





1 


32.50 


1 




0.045 53 


8.000 


0.2222 


- 


K M 


MP19 


1 


0,±1 


-1 


-47.58 


1 




0.076 13 


-9.795 


0.2222 


~ 


K M 


MP20 


1 


1 





1 


31.77 


0.016 75 




1 


0.2832 


~ 


M K 


MP21 


1 


1 





1 


24.98 


0.010 90 




1 


0.2825 


~ 


M K 


MP22 


1 


-1 


0,±1 


1 


-0.8418 


0.8922 




1 


0.4325 


~ 


M K 


MP23 


1 


±1 


-1 





1 




1 


0.7778 


0.2222 


~ 


S± M 


MP24 


1 


-1 


±1 


1 





1 




1 





~ 


M S± 


MP25 


1 


1 


1 


1 


1 


0.2210 


0.3839 


1 


0.2222 


~ 


c c 


MP26 


1 


-1 


1 


1 


1 


0.069 37 


0.3495 


1 


0.2222 


~ 


M C 


MP27 


1 


-1 


-1 


1 


1 


1.296 


0.6584 


1 


0.2222 


~ 


M M 


MP28 


1 


-1 


-1 


1 


1 


0.5507 


0.8095 


1 


0.2222 


~ 


M M 



Table 11: Exact solutions for 5 = —0.1. e = and 1 correspond to generalized de Sitter solutions (a ~ e'^*, & ~ e"') 
and power law ones (a ~ r'', b ^ r"), respectively. A is a power exponent of power law solutions in the Einstein frame 
{as - 4)-K, S±, So and M mean a kinetic dominance, a static space with positive (or negative) curvature, a flat static 
space, and a Milne-type space, respectively. ME3± and ME6+ are inflationary solution in the Einstein frame, but our 
3-space is shrinking in ME3_ . 



Solution 


e 


CTp 


(Tg 




1/ 


Oo 


bo 


A 




Type 


ME3± 











±0.1682 


±0.1682 






1.286 


±0.2857 


K K 


ME4± 











±0.7429 


TO.4540 






0.5326 


±0.2857 


K K 


ME6± 








1 


±0.3072 







4.605 







K S+ 


ME9± 





1 








±0.2011 


2.658 




1 


±0.2857 


S+ K 


ME10± 





-1 








±0.3141 


1.149 




1 


±0.2857 


S- K 


ME12 


1 





-1 





1 




1 


0.7778 


0.2222 


So M 


ME13 


1 


-1 





1 





1 




1 





M So 
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